Renewable  and  Sustainable  Energy  Reviews  15(2011)  1 893-1917 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Renewable  and  Sustainable  Energy  Reviews 

journal  homepage:  www.elsev er.com/locate/rser 


Mathematical  modeling  of  solid  oxide  fuel  cells:  A  review 

S.  Ahmad  Hajimolana* 1 2 3,  M.  Azlan  Hussain3  *,  W.M.  Ashri  Wan  Daud3, 

M.  Soroush  b,  A.  Shamiri3 

a  Department  of  Chemical  Engineering,  University  of  Malaya,  50603  Kuala  Lumpur,  Malaysia 
b  Department  of  Chemical  and  Biological  Engineering,  Drexel  University,  Philadelphia,  PA  19104,  USA 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  23  November  2010 
Accepted  27  December  2010 


Keywords: 

Solid  oxide  fuel  cell 

Mathematical  modeling 

Dynamics 

Heat  transfer 

Mass  transfer 

Electrochemistry 

Polarization 

Novel  SOFC  models 


This  paper  presents  a  review  of  studies  on  mathematical  modeling  of  solid  oxide  fuel  cells  (SOFCs)  with 
respect  to  the  tubular  and  planar  configurations.  In  this  work,  both  configurations  are  divided  into  five 
subsystems  and  the  factors  such  as  mass/energy/momentum  transfer,  diffusion  through  porous  media, 
electrochemical  reactions  with  and  without  CO  oxidation,  shift  and  reforming  reactions,  and  polarization 
losses  inside  the  subsystems  are  discussed.  Using  variety  of  fuels  fed  to  SOFCs  is  issued  and  their  effect 
on  the  system  is  compared  briefly.  A  short  review  of  solid  oxide  fuel  cell  configurations  and  different 
flow  manifolding  are  also  presented  in  this  study.  Novel  models  based  on  statistical  data-driven  approach 
existing  in  the  literatures  are  considered  shortly.  Although  many  studies  on  solid  oxide  fuel  cells  modeling 
have  been  done,  still  more  research  needs  to  be  done  to  improve  the  models  in  order  to  predict  the  fuel  cell 
behaviors  more  accurately.  At  the  end  of  this  paper  the  works  and  studies  that  can  be  done  for  improving 
the  fuel  cell  models  is  suggested  and  pointed  by  the  authors. 
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1.  Introduction 

Ever  increasing  energy  consumption,  rising  public  awareness 
for  environmental  protection  and  higher  prices  of  fossil  fuels  have 
motivated  many  to  look  for  alternative/renewable  energy  sources. 
Fossil  fuel  resources  are  limited  and  are  expected  to  end  early  in 
the  next  century.  By  2015,  the  world  fossil  fluid  fuel  demand  will 
exceed  the  world  fossil  fluid  fuel  production,  which  is  expected  to 
lead  to  an  energy  shortage  crisis  unless  a  sustainable  alternative 
fuel  will  be  available  by  them  [1-3]. 

Small-scale  power-generation  systems  such  as  wind  turbines, 
photovoltaics,  micro-turbines  and  fuel  cells  can  play  an  important 
role  in  meeting  consumer  demands  for  “greener”  energy.  These 
small-scale  generation  systems,  usually  located  near  customers 
rather  than  at  central  or  remote  locations,  can  be  used  to  form 
distributed  power  generation  networks.  Among  different  types  of 
small-scale  power  generation  systems,  fuel  cells  have  received 
more  attention,  because  they  can  provide  both  heat  and  power. 
Household  fuel  cell  lets  households  to  generate  part  of  the  electric¬ 
ity  they  require,  reducing  the  electricity  purchase  from  an  electric 
power  company.  In  addition,  they  serve  as  a  cogeneration  system, 
which  use  thermal  energy  produced  by  electricity  generation  for 
water  and  air  heating.  It  would  be  best  if  exhaust  heat  emitted  dur¬ 
ing  the  generation  process  from  large-scale  power  stations  could 
also  be  converted  to  thermal  energy  for  heating  air  or  water  to  use 
the  inputted  energy  more  efficiently.  However,  it  is  not  considered 
practical  to  send  hot  water  (thermal  energy)  to  homes  by  pipelines 
from  far-removed  locations  of  such  large-scale  power  stations. 

Fuel  cells  are  static  energy  conversion  devices  that  partially  con¬ 
vert  the  chemical  energy  of  fuels  directly  into  electrical  energy  and 
produces  water  as  its  byproduct  [4,5].  The  basic  principle  of  the  fuel 
cell  was  discovered  in  the  year  1838  by  Swiss  scientist  Christian 
Friedrich  Schnbein.  In  1839  Sir  William  Grove  developed  the  first 
fuel  cell  based  on  reversing  the  electrolysis  of  water  by  an  accident 
[3].  In  1950  Francis  Bacon  at  Cambridge  University  demonstrated 
the  first  5  kW  alkalinahoe  fuel  cell.  After  the  successful  develop¬ 
ment  of  alkaline  fuel  cells,  NASA  needed  a  compact  system  to 
generate  electricity  for  space  shuttle  applications.  In  1970s,  inter¬ 
national  fuel  cells  developed  a  12kW  alkaline  fuel  cell  for  NASA’s 
space  shuttle  orbiter  to  supply  reliable  power  without  the  use  of 
any  backup  powers  like  batteries.  Beginning  in  the  mid-1960s,  the 
research  work  was  focused  on  further  development  of  various  fuel 
cells  for  applications  like  stationary  powers  and  transportations. 
Further  the  government  agencies  in  the  USA,  Canada,  Japan  and 
Malaysia  have  significantly  increased  their  funding  for  fuel  cell  in 
R&D  [6]. 

There  are  six  main  types  of  fuel  cell:  alkaline  fuel  cell  (AFC), 
phosphoric  acid  fuel  cell  (PAFC),  molten  carbonate  fuel  cell 
(MCFC),  polymer  electrolyte  membrane  fuel  cell  (PEMFC)  and  direct 
methanol  fuel  cell  (DMFC),  solid  oxide  fuel  cell  (SOFC).  Table  1 
depicts  their  characteristics  briefly. 

After  the  development  of  power  conversion  devices,  much  more 
researches  are  doing  on  solid  oxide  fuel  cells  to  reduce  its  higher 
installation  cost.  Solid  oxide  fuel  cells  are  considered  to  be  one  of  the 
most  advanced  designs  for  mid-  to  large-scale  applications.  They 
are  among  the  most  promising  types  of  fuel  cells  currently  being 
considered  as  a  power  source  for  automobiles  and  stationary  power 


plants  [2,3,7].  Because  the  electrolyte  is  a  layer  of  ceramic  material 
with  high-temperature  endurable  porous-media  electrodes,  SOFCs 
can  generally  operate  in  a  high-temperature  range  (800-1000  C). 
High  operating  temperature  has  some  advantages,  such  as  high 
energy  conversion  efficiency,  flexibility  of  usable  fuel  type,  and  high 
temperature  exhaust  gas.  Disadvantages  include  potential  ther¬ 
mal  fatigue  failure  of  the  cell  material  and  sealing  under  the  high 
temperature,  as  well  as  the  fact  that  cell  temperature  fluctuations 
induce  thermal  stress  in  the  cell  ceramics.  Thus,  it  is  important  to 
operate  SOFCs  in  such  a  way  that  the  stack  temperature  remains 
within  a  tight  design  range. 

In  this  case,  modeling  is  a  low  cost  method  for  studying  and 
investigating  the  fuel  cell  in  order  to  optimize  and  control  the  fuel 
cell  behavior,  increase  its  efficiency  and  performance,  and  finally 
decrease  the  higher  installation  cost.  Among  many  mathematical 
modeling  studies,  there  are  some  review  papers,  concerning  differ¬ 
ent  aspects  of  solid  oxide  fuel  cells  investigation  which  introduce  a 
complete  view  of  all  the  proceeding  works  to  the  researchers  who 
are  going  to  develop  and  improve  them.  In  this  regard,  Badwal  et  al. 
[8]  reviewed  the  different  kinds  of  materials  used  in  electrolyte  and 
electrodes.  Similar  work  has  been  done  by  Zhu  et  al.  [9],  presented  a 
review  paper  about  anode  materials  used  in  SOFCs.  The  criteria  for 
the  anode  of  SOFC  were  first  discussed.  They  also  discussed  about 
the  problems  existing  in  the  developed  anode  materials.  In  this 
paper  the  advantages  and  disadvantages  of  anodes  materials  were 
illustrated  in  order  to  give  a  view  for  further  research  and  develop¬ 
ment  of  new  generation  of  anode  materials  for  SOFC.  A.  Stambouli 
[6]  reconsidered  SOFCs  as  environmentally  clean  and  an  efficient 
source  of  energy.  The  advantages  and  disadvantages  of  models 
regarding  to  AC  impedance  was  studied  completely  with  Q.  Huang 
et  al.  [10].  Zhang  et  al.  [1 1  ]  performed  a  literature  study  about  dif¬ 
ferent  strategies  and  concepts  for  SOFC-based  integration  systems, 
which  are  timely  transformational  energy-related  technologies 
available  to  overcome  the  threats  posed  by  climate  change  and 
energy  security.  A  complete  review  paper  about  dynamic  models  of 
single-cell  SOFCs  with  respect  to  various  methods,  assumptions  and 
models  in  design,  diagnosis  and  SOFC  operation  was  presented  by 
Bhattacharyya  et  al.  [12]  but  it  was  less  discussion  on  steady-state 
models.  The  impact  of  the  cell  geometry  was  also  presented  in  their 
review  paper.  S.  Kaka^  et  al.  [  1 3  ]  reconsidered  the  numerical  models 
of  SOFCs  and  summarized  the  present  status  of  the  SOFC  modeling. 
V.  Lawlor  et  al.  [14]  discussed  in  two  parts  of  a  review  paper  about 
micro-tubular  (MT-SOFC)  solid  oxide  fuel  cell.  Part  I  illustrated  the 
reader  to  the  MT-SOFC  stack  and  its  applications,  indicating  who  is 
researching  what  in  this  field  and  also  specifically  investigate  the 
design  issues  related  to  multi-cell  reactor  systems  called  stacks. 
Part  II  [15]  reviewed  in  detail  the  combinations  of  materials  and 
methods  used  to  produce  the  electrodes  and  electrolytes  of  MT- 
SOFC’s.  Also  the  role  of  modeling  and  validation  techniques  used 
in  the  design  and  improvement  of  the  electrodes  and  electrolytes 
was  investigated.  M.  Andersson  et  al.  [16]  presented  a  literature 
study  with  respect  to  various  transport  processes  such  as  mass, 
heat,  momentum,  and  also  charge,  electrochemical,  shift  and  inter¬ 
nal  reforming  reactions.  This  work  focused  on  SOFCs  designing  by 
combining  the  accuracy  at  microscale  with  the  calculation  speed  at 
macroscale,  based  on  clear  understanding  of  transport  phenomena, 


Table  1 

Types  of  fuel  cell  [1  ]. 
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Fuel  cell 

Temperature  (°C) 

Efficiency  (%) 

Application 

Advantage 

Disadvantages 

Alkaline  fuel  cell  (AFC) 

50-90 

50-70 

Space  application 

High  efficiency 

Intolerant  to  CO2  in  impure  H2 
and  air,  corrosion,  expensive 

Phosphoric  acid  fuel  cell  (PAFC) 

175-220 

40-45 

Stand-alone  &  combined 
heat  &  power 

Tolerant  to  impure  H2 

Low  power  density,  corrosion 
&  sulfur  poisoning 

Molten  carbonate  fuel  cell  (MCFC) 

600-650 

50-60 

Central,  stand-alone  & 
combined  heat  &  power 

High  efficiency 

Electrolyte  instability, 
corrosion  &  sulfur  poisoning 

Solid  oxide  fuel  cell  (SOFC) 

800-1000 

50-60 

Central,  stand-alone  & 
combined  heat  &  power 

High  efficiency  &  direct 
fossil  fuel 

High  temperature,  thermal 
stress  failure,  coking  &  sulfur 
poisoning 

Polymer  electrolyte  membrane 
fuel  cell  (PEMFC) 

60-100 

40-50 

Vehicle 

High  power  density,  low 
temperature 

Intolerant  to  CO  in  impure  H2 
and  expensive 

Direct  methanol  fuel  cell  (DMFC) 

50-120 

25-40 

No  reforming,  high  power 
density  &  low  temperature 

Low  efficiency,  methanol 
crossover  &  poisonous 
byproduct 

chemical  reactions  and  functional  requirements.  Useful  methods 
were  investigated  for  SOFCs  modeling.  Coupling  methods  between 
different  approaches  and  length  scales  by  multiscale  models  were 
outlined  in  order  to  enhance  the  understanding  for  detailed  trans¬ 
port  phenomena  and  making  a  correct  decision  on  the  specific 
design  and  control  of  operating  conditions.  It  was  remarked  that 
the  connection  between  numerical  modeling  and  experiments  is 
too  rare  and  also  that  material  parameters  in  most  cases  are  valid 
only  for  standard  materials  and  not  for  the  actual  SOFC  component 
microstructures. 

A  review  study  on  mathematical  modeling  based  on  steady- 
state  and  dynamic  behavior,  and  control  of  polyelectrolyte 
membrane  and  SOFCs  with  respect  to  zero-  one-,  two-,  and  three- 
dimensional  models  was  performed  by  M.  Bavarian  et  al.  [17].  The 
behavior  of  the  fuel  cell  and  processes  that  occur  inside  the  fuel 
cells  and  contribute  the  existence  of  multiple  time-scales  in  the  fuel 
cells  were  investigated  and  also  important  components  of  these 
models  were  remarked.  It  was  discussed  that  under  which  situa¬ 
tions  a  fuel  cell  observes  steady  state  multiplicity  and  also  stability 
of  the  steady  state  was  issued.  In  continue,  a  review  of  control 
configurations  and  strategies  by  mentioning  the  advantages  and 
disadvantages  of  each  was  presented.  Finally,  the  topics  that  require 
further  research  studies  are  discussed. 

In  this  study,  mathematical  modeling  of  solid  oxide  fuel  cells 
with  respect  to  the  tubular  and  planar  configurations  is  reviewed 
into  five  subsystems  considering  the  factors  such  as  polarization 
losses,  mass/energy/momentum  conservations,  diffusion  through 
porous  media,  electrochemical  phenomena  in  the  PEN  region  and 
shift/reforming  reactions  inside  them  which  is  novel  in  this  regard. 
Also,  using  variety  of  fuels  fed  to  the  SOFCs  is  discussed  and  their 
effect  on  the  system  is  compared  briefly.  A  short  review  of  solid 
oxide  fuel  cell  configurations  and  different  flow  manifolding  are 
also  presented  in  this  study.  For  the  first  time  in  this  area,  novel 
models  based  on  statistical  data-driven  approach  existing  in  the 
literatures  are  studied.  At  the  end  of  the  paper,  it  is  suggested  and 
pointed  which  topics  need  more  study  for  improving  the  fuel  cell 
models. 

The  remaining  part  of  this  article  is  organized  as  follows:  Sec¬ 
tion  2  presents  general  principle  of  solid  oxide  fuel  cell;  Section  3 
discusses  mathematical  modeling  of  SOFCs,  Section  4  addresses  lit¬ 
erature  on  the  SOFC  configuration;  and  Conclusion  is  presented  in 
Section  5. 


2.  General  principle  of  solid  oxide  fuel  cell 

A  solid  oxide  fuel  cell  uses  a  hard  ceramic  electrolyte  and  oper¬ 
ates  at  temperatures  up  to  1000°C.  A  mixture  of  zirconium  oxide 
and  calcium  oxide  forms  a  crystal  lattice,  although  other  oxide  com¬ 


binations  have  also  been  used  as  electrolytes.  The  crystal  lattice 
allows  oxide  ions  to  pass  through  it  to  reach  the  anode  surface, 
where  the  oxide  ions  combine  with  H+  ions  and  form  water.  The 
solid  electrolyte  is  coated  on  both  sides  with  specialized  porous 
electrode  materials.  The  anode  consists  of  a  metallic  nickel-  and 
Y203-stabilized  ZrC>2  skeleton,  which  inhibits  sintering  of  the  metal 
particles  and  provides  a  thermal  expansion  coefficient  compara¬ 
ble  to  those  of  the  other  cell  materials,  thus  limiting  the  buildup 
of  stresses  resulting  from  a  difference  in  the  coefficient  of  ther¬ 
mal  expansion.  The  anode  structure  has  a  porosity  of  20-40%  to 
facilitate  mass  transport  of  the  reactant  and  product  gases.  The 
cathode  material  mostly  used  is  strontium-doped  lanthanum  mag¬ 
netite  (Lai_xSrxMn03,  for  x  0.10-0.15);  its  structure,  such  as  that  of 
the  anode,  is  porous  to  permit  rapid  mass  transport  of  the  reac¬ 
tant  and  product  gases.  The  cathode  material  has  low  levels  of 
chemical  reactivity  with  the  electrolyte,  which  extends  the  life¬ 
time  of  the  material.  However,  it  is  a  poor  ionic  conductor,  so 
the  electrochemically  active  reactions  are  limited  to  the  triple¬ 
phase  boundary  (TPB),  where  the  electrolyte,  gas,  and  electrode 
meet.  In  other  words,  the  TPB  of  a  fuel  cell  is  the  area  of  contact 
between  the  three  phases  necessary  for  electrochemical  (hydro¬ 
gen  oxidation  and  oxygen  reduction)  reactions  at  the  electrode: 
an  ion-conducting  phase,  an  electron  conducting  phase,  and  a  gas 
phase  [18].  The  bigger  the  TPB  area  of  a  cell,  the  better  the  qual¬ 
ity  of  the  cell;  a  bigger  TPB  area  allows  the  reactions  to  occur 
in  more  sites,  thus  maximizing  current  flow.  Lanthanum  stron¬ 
tium  magnetite  also  works  well  as  a  cathode  at  high  temperatures, 
but  its  performance  quickly  decreases  as  operating  temperature 
decreases  to  less  than  80  °C.  The  electrochemical  reaction  occurs 
when  both  fuel  and  air  are  allowed  to  flow  into  the  SOFC  fuel  cell 
separately.  They  diffuse  through  the  porous  electrode  structure  to 
the  interlayer  and  are  adsorbed.  At  the  cathode  interlayer,  oxygen 
is  reduced  by  the  incoming  electrons  to  produce  oxygen  anions 
that  are  conducted  through  the  electrolyte  to  the  anode  inter¬ 
layer  where  they  electrochemically  combine  with  the  adsorbed 
hydrogen  to  form  water  and  release  electrons  to  the  external  cir¬ 
cuit. 

The  main  advantages  of  the  SOFC  is  that  they  can  be  operated  at 
high  efficiency  of  50-60%  and  a  separate  reformer  is  not  required 
to  extract  hydrogen  from  the  fuel  due  to  its  internal  reforming 
capability.  Waste  heat  can  be  recycled  to  make  additional  elec¬ 
tricity  by  cogeneration  operation  [19,20].  The  slow  start  up,  high 
cost  and  intolerant  to  sulfur  content  of  the  fuel  cell  are  some  of  its 
drawbacks.  It  is  not  suitable  for  larger  fluctuations  in  load  demand. 
Therefore,  the  SOFC  is  mainly  used  for  medium  and  large  power 
applications.  In  1997  the  Ceramic  Fuel  Cells  Limited  Company  was 
demonstrated  a  5  kW  laboratory  prototype  fuel  cell  system.  Yakabe 
et  al.  [21]  developed  a  3  kW  SOFC  at  Tokyo  gas  Co.  Ltd.  And  they 
also  analyzed  the  key  factors  to  improve  the  performance  of  SOFC 
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Table  2 

The  advantages  and  disadvantages  of  PSOFC  and  TSOFC. 


Type  of  SOFC 

Advantages 

Disadvantages 

PSOFC 

-Has  the  potential  to  offer  higher  power  density  than  the  tubular 
design. 

-  Due  to  its  compactness,  it  can  be  stacked  in  resemblance  to  polymer 
electrolyte  membrane  (PEM)  fuel  cells  to  satisfy  the  power 
requirement  of  an  application. 

-  Planar  SOFC  is  simple  to  fabricate  and  can  be  manufactured  into 
various  configurations. 

-The  problems  associated  with  the  development  and 
commercialization  of  planar  SOFCs  is  requirement  of  high 
temperature  gas  seals. 

-  Internal  stresses  in  cell  components  due  to  non-uniform 
temperature  distributions  and  high  manufacturing  cost. 

TSOFC 

-  Show  great  promise  as  potential  candidates  to  replace  conventional 
heat  engines. 

-Have  higher  energy  efficiency. 

-  They  are  designed  for  mid-  to  large-scale  applications  up  to  2  MW. 

-Lower  power  density  than  the  PSOFC  design. 

-  Have  high  ohmic  loss 

in  the  micro-grid  system  and  are  building  a  250  kW  commercial 
stack  model. 

2.1.  SOFC  configurations 

There  are  two  principal  SOFC  configurations,  tubular  and  pla¬ 
nar.  The  Schematic  of  tubular  and  planar  SOFC  is  illustrated  in 
Figs.  1  and  2  respectively.  Tubular  configuration  consists  of  two 
tubes,  an  outer  tube  and  an  inner  tube.  The  outer  tube  is  the  cell 
tube.  The  outer  surface  of  the  cell  tube  is  anode  side  of  the  cell, 
and  its  inner  surface  is  the  cathode  side.  Between  the  anode  and 
cathode  sides  lies  the  solid  oxide  electrolyte.  The  inner  tube  is  an 
air  injection  and  guidance  tube,  composed  of  alumina,  from  which 
preheated  air  is  injected  into  the  bottom  of  cell  tube  and  flows  over 
the  cathode  surface  of  the  cell  tube  through  the  gap  between  the 
injection  and  cell  tubes.  The  end  of  the  cell  tube  is  closed.  Fuel  gas 
flows  over  the  anode  surface  through  the  gap  among  the  cell  tubes. 
Oxygen  ions  pass  through  the  cathode  and  electrolyte  and  react 
with  fuel,  creating  an  electric  current.  A  major  advantage  of  tubu¬ 
lar  Solid  Oxide  Fuel  Cells  (TSOFCs)  over  other  types  of  fuel  cells  is 
that  a  variety  of  hydrocarbon-based  gases  or  their  synthesis  deriva¬ 
tives,  such  as  natural  gas,  biomass  and  coal  can  be  potentially  used 
as  fuel  sources.  CO  present  in  the  synthesis  gas  can  be  oxidized  in 
TSOFCs  to  generate  C02  and  electrical  energy.  The  disadvantage  of 
tubular  SOFC  compared  to  planar  SOFC  is  that  the  ohmic  losses  are 
much  higher  due  to  the  longer  ways  of  electrical  current. 

Generally,  planar  configuration  is  constructed  of  a  positive- 
electrolyte-negative  electrode  (PEN),  interconnect  top  and  bottom 
of  the  cell  stack,  air  and  fuel  channels.  The  fuel  gas  channel  is 
between  the  anode  and  the  separator  plate  and  on  the  other  side 
the  air  channel  is  located  between  the  cathode  and  separator  plate. 
Recognizing  the  advantages  of  the  low-cost,  simplicity  in  manu¬ 
facturing  and  high  volume  manufacturing  with  high  volumetric 
power  densities,  most  fuel  cells  manufacturers  are  concentrating 
on  the  planar  SOFC  concept.  However,  a  challenge  with  the  planar 
geometries  is  in  obtaining  mechanically  stable  structure,  as  thin 
layer  ceramics  are  inherently  susceptible  to  failure  when  subjected 
to  moderate  stresses.  The  planar  concept  provides  the  largest  area 
for  the  current  collector  but  it  shows  disadvantages  in  start-up  time 
and  mechanical  stability.  The  advantages  and  disadvantages  of  two 
types  of  SOFC  are  depicted  in  Table  2. 

Based  on  the  thicknesses  of  the  electrodes  and  electrolyte 
layers,  SOFCs  are  divided  also  into  three  types  namely  anode- 
supported,  electrolyte-supported  and  cathode-supported.  It  was 
explored  that  ohmic  and  polarization  losses  in  electrode-supported 
are  smaller  than  in  electrolyte-supported  cells  [22]. However,  the 
anodic  polarization  is  smaller  than  the  cathodic  one  [23].  More¬ 
over,  anode-supported  solid-oxide  fuel-cells  can  operate  at  lower 
temperatures,  which  results  in  diminished  manufacturing  costs 
[24].  For  these  reason  anode-supported  cells  are  often  considered 
to  be  the  preferred  configuration  for  SOFCs.  But  the  thick  anode 
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Fig.  1.  Schematic  of  tubular  SOFC. 


supported  structure  leads  to  enhance  the  concentration  polariza¬ 
tion  due  to  diffusion  limitation,  and  thus  decreases  the  fuel  cell 
performance  [25]. 

Regarding  the  flow  configuration,  there  are  three  main  types 
flow  inside  the  fuel  cells  namely  co-flow,  counter-flow  and  cross- 
flow. 


3.  Mathematical  modeling  of  solid  oxide  fuel  cells 

Mathematical  modeling  is  an  essential  tool  in  designing  fuel  cell 
systems  which  represents  the  important  aspects  of  the  existing 
system  and  presents  knowledge  of  that  system  in  usable  form.  In 
other  word,  such  models  can  provide  a  picture  of  voltage,  current 
density,  temperature,  velocity  and  concentration  of  components 
as  functions  of  position  and  time  for  various  cell  configurations 
and  operating  conditions.  The  effect  of  different  variables  and  their 
sensitivity  on  the  fuel  cell  performance  can  be  studied  by  these 
mathematical  models.  A  model  is  also  useful  in  predicting  the 
effects  of  altering  process  variables  and  in  using  that  information 
to  optimize  cell  performance.  Moreover,  dynamic  simulation  can 
provide  reference  for  the  control  system  research  and  design  in 
future. 

Solid  oxide  fuel  cell  modeling  for  both  the  tubular  and  planar 
configurations  tacked  into  account  the  electrochemical  model,  dif¬ 
fusion  through  porous  media,  the  mass,  energy,  and  momentum 
conservation  balances. 

To  perform  first-principles  model  of  the  SOFC  for  the  both  con¬ 
figurations,  single  tubular  and  planar  fuel  cells  are  considered  and 
divided  into  several  subsystems  as  shown  in  Figs.  3  and  4. 

For  the  tubular  SOFC  the  subsystems  are  as  follow: 

•  Subsystem  1  (SSI):  air  inside  the  injection  tube; 

•  Subsystem  2  (SS2):  injection  tube; 
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Fig.  2.  Schematic  of  planar  SOFC. 


Although  mass  transfer  via  diffusion  inside  the  fuel  cell  sig¬ 
nificantly  affects  the  fuel  cell  performance,  some  researchers  like 
Costamagna  et  al.  [26],  Burt  et  al.  [27],  Stiller  et  al.,  [28],  Ji  et  al.  [29] 
neglected  diffusion  of  components  through  media  electrodes.  Burt 
et  al.  [27]  were  assumed  that  changes  in  the  x-direction  are  small 
and  can  safety  be  ignored. 

The  diffusion  through  the  porous  media  was  modeled  by  writing 
detailed  mass  transfer  equations  with  many  researchers  like  Chan 
et  al.  [30],  Yakabe  et  al.  [31  ],  Bove  et  al.  [32],  Ho  et  al.  [33],  Qi  et  al. 
[34]  and  others  [35-41].  In  these  works,  it  was  observed  that  the 
fuel  cell  behavior  was  more  accurately  than  the  models  neglected 
the  diffusion  of  gases  in  the  electrodes. 

The  diffusion  of  component  through  porous  electrodes  is 
assumed  by  using  one  of  the  three  approaches  as  the  following: 


•  Subsystem  3  (SS3):  air  inside  the  space  between  the  cell  and 
injection  tubes; 

•  Subsystem  4  (SS4):  cell  tube; 

•  Subsystem  5  (SS5):  fuel  mixture  over  the  anode  side  of  the  cell 
tube. 


(a)  Fields  Law  model  (FLM)  is  the  simplest  diffusion  model  and  used 
in  dilute  or  binary  systems  by  Stiller  et  al.  [28],  Chan  et  al.  [30], 
Bove  et  al.  [32],  Hussain  et  al.  [37]  and  many  other  researchers 
[30,32,37-40,42-53].  The  model  is  given  by: 

Ni  =  [DfV(pYf)]  (4) 


For  the  planar  SOFC  the  subsystems  are  as  follow: 

•  Subsystem  1  (SSI):  air  channel; 

•  Subsystem  2  (SS2):  air  side  interconnector; 

•  Subsystem  3  (SS3):  fuel  side  interconnector; 

•  Subsystem  4  (SS4):  cell  tube; 

•  Subsystem  5  (SS5):  fuel  channel. 

The  SOFCs  models  are  then  derived  by  writing  the  appropriate 
electrochemical  and  transfer  phenomena  equations. 


3.2.  Model  equations  for  subsystem  2  (SSI) 

This  subsystem  for  TSOFC  includes  the  air  inside  the  injec¬ 
tion  tube  and  for  PSOFC  is  the  air  flow  inside  the  air  channel, 
between  air  side  interconnector  and  cathode.  In  this  subsystem, 
mass/energy/momentum  transfer  models  can  be  accounted. 

3.2.2.  Mass  transfer  model 

Gas  transport  is  important  to  enhance  the  performance  of  the 
fuel  cell.  The  performance  of  a  SOFC  at  high  operating  current 
densities  is  reduced  mainly  due  to  the  mass  transport  losses. 
Micro-structural  parameters  also  largely  affect  the  gas  transport. 
Mass  transport  models  inside  the  electrodes  are  giving  an  accurate 
prediction  for  gas  concentration  at  the  anode-electrolyte-cathode 
interface.  Therefore;  an  accurate  modeling  of  this  phenomenon  is 
imperative  for  development  of  better  fuel  cell  designs. 


(i)  A  mass  balance  on  the  air  inside  SSI  for  TSOFC  takes  the  form 


dpin)  .... 

+  v(umv2)  = 0 


(i) 


(ii)  A  mass  balance  on  the  air  and  02  inside  SSI  (cathode  channel) 
for  PSOFC  takes  as  the  follow: 


dncKc 

-yjr + v(u«?;Va?;c)  =  -™0; 


dp 


ch,c 

air 


dt 


Y°l+V(uch’cpch’cY0,)  = 

y  air  rair  u2 J 


VNo- 


(2) 


(3) 


where  i  denotes  the  generic  ith  specie,  y02  the  mass  fraction  of 
the  02  species,  p  is  the  density,  u  is  the  gas  velocity  and  VNq2 
is  the  diffusion  of  02  into  the  porous  cathode  for  participat¬ 
ing  in  electrochemical  reaction,  Superscripts  inj  and  ch,c  denote 
injection  tube  and  cathode  channel,  respectively. 


where  D,  is  the  diffusion  coefficient. 

(b)  The  Stefan-Maxwell  model  (SMM)  is  more  commonly  used  in 
multi-component  systems  for  calculation  of  specie  fluxes.  The 
model  is  given  by: 


n 

-V(Y,.)=  £ 


ytNi  -  YfNj 

-ai— 


where  N  is  the  species  flux  and  C  is  the  concentration. 

Some  other  researchers  like  Bhattacharyya  et  al.  [54],  Virkar 
et  al.  [55]  Kim  et  al.  [56]  Bove  et  al.  [57],  Sanchez  et  al.  [58,59] 
and  Jin  et  al.  [60]  considered  SMM  model  in  their  survey. 

(c)  The  Dusty  Gas  model  (DGM),  is  the  extended  model  of  Stefan- 
Maxwell  equation,  which  is  the  most  complete  diffusion  model 
considering  molecule-molecule  and  molecule-wall  interactions 
in  the  porous  media,  commonly  used  as  well,  because  it  takes 
into  account  Knudsen  diffusion  phenomena  and  neglect  the 
assumption  of  equal-molar  counter  diffusion  in  Fields  law, 
which  becomes  invalid  if  the  molecular  weights  of  the  diffusing 
gas  species  differ  widely.  The  model  is  given  by: 


Nf 


n 


D 


+ 


i,k 


E 

j=i,j  ~t~  i 


m 


YiNj 


CDi 


hj 


In  many  surveys  this  model  is  used  [38,49,61-71].  Tseronis 
et  al.  [72]  and  Cayan  et  al.  [40]  used  a  combination  of  diffusion 
method. 

Tseronis  et  al.  [72]  discussed  multidimensional  model  trans¬ 
port  phenomena  in  an  anode  channel  and  diffusion  through  porous 
anode  in  an  isothermal  planer  SOFC  based  on  a  combination  of  the 
Stefan-Maxwell  and  Dusty  gas  method.  The  results  were  validated 
with  experimental  data  (Yakabe  et  al.  [31  ]). 

Cayan  et  al.  [40]  illustrated  a  two-dimensional  model  based  on 
mass  transfer  inside  an  SOFC.  In  this  survey,  diffusion  methods, 
Stefan-Maxwell  and  Fields  law  relations,  both  including  Knudsen 
diffusion,  in  order  to  predict  the  species  concentration  gradient  was 
compared.  The  results  showed  that  at  low  current  densities,  both 
models  are  well  in  agreement  while  by  increasing  the  current  the 
differences  becomes  bigger. 

A  complete  survey  for  performance  comparison  of  Fields  model, 
dusty  gas  model  and  Stephan-Maxwell  model  have  been  done  by 
R.  Suwanwarangkul  et  al.  [52].  It  was  showed  that  DGM  is  the 
most  suitable  model  for  the  H2-H20  and  CO-C02  systems  due  to 
the  Knudson  diffusion  impact.  However,  due  to  its  complexity,  this 
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model  can  only  be  used  at  the  conditions  of  low  reactant  concen¬ 
tration,  high  operating  current  density  and  small  pore  size  where 
high  accuracy  of  model  prediction  is  needed.  Otherwise,  the  SMM 
is  a  useful  model  for  H2-H20  system  and  FM  is  suitable  for  the 
C0-C02  system.  SMM  does  not  tacked  into  account  the  Knudson 
diffusion  term  in  order  to  calculate  the  pore  size  effect  while  both 
FM  and  DGM  do.  It  was  fund  that  the  equimolar  counter  diffusion 
is  assumed  to  determine  the  flux  ratio  in  both  FM  and  the  SMM.  By 
contrast,  the  flux  ratio  in  the  DGM  is  calculated  from  the  ratio  of 
the  square-root  of  the  gas  molecular  weight.  It  was  illustrated  that 
DGM  is  suitable  for  multi-component  system  of  H2-H20-C0-C02, 
because  the  flux  ratio  calculation  is  complicated  and  equimolar 
counter  diffusion  cannot  be  assumed.  Moreover,  the  Knudson  dif¬ 
fusion  effect  must  be  considered. 

Diffusion  through  the  porous  material  is  typically  described  by 
either  ordinary  or  Knudsen  diffusion  and  has  been  found  to  play 
an  important  role  in  catalytic  reaction.  Ordinary  diffusion  occurs 
when  the  pore  diameter  of  the  material  is  large  in  comparison 
to  the  mean  free  path  of  the  gas  molecules.  Molecular  transport 
through  pores  which  are  small  in  comparison  to  the  mean  free 
path  of  the  gas  is  regarded  as  a  Knudsen  type  of  diffusion.  For 
Knudsen  diffusion,  molecules  collide  more  frequently  with  the  pore 
walls  than  with  other  molecules.  Upon  collision,  the  atoms  are 
instantly  adsorbed  on  to  the  surface  and  are  then  desorbed  in  a 
diffusive  manner.  As  a  result  of  frequent  collisions  with  the  wall 
of  the  pore,  the  transport  of  the  molecule  is  impeded.  The  Knud¬ 
sen  diffusion  coefficient  can  be  predicted  using  kinetic  theory  by 
relating  the  diameter  of  the  pore  and  the  mean  free  path  of  the  gas. 
For  straight  and  round  pores,  the  diffusion  coefficient  is  given  by 
[30]: 


97.01* 


where  Tct  is  the  cell  tube  temperature,  and  r  is  mean  pore  radius 
which  can  be  calculated  from  the  bulk  density,  SA  area  of  the 
porous  solid,  and  s  porosity  [30]. 


2s 

SaPb 


To  account  for  the  tortuous  path  of  the  molecule  rather 
than  along  the  radial  direction  and  the  porosity  of  the  solid 
electrode  for  the  fact  that  diffusion  occurs  only  in  the  pore 
space,  an  effective  Knudsen  diffusion  coefficient  may  be  used,  i.e. 
[30]: 


Neufield  et  al.  [73]  is  given  by: 

1.06036  0.19300 

1  "  (kr/e^)015610  +  exp(0.47635(kr/e/,B)) 

1.03587  1.76474 

+  exp(1.52996(kT/eAB))  +  exp(3.8941  l(kr/eAB)) 


(12) 


Hajimolana  and  Soroush  [22],  Qi  et  al.  [34],  Bhattacharyya  et  al. 
[54]  and  other  authors  [22,34,53,54,59,64,77-79]  used  Fuller  et  al. 
correlation  which  is  given  by  [80]: 


1.013  x  10-7rcV75[(l/Mo2)  +  (l/MN2)]1/2 

K2B + pXo  -  s0a()  m/sAB)2 


(13) 


The  effective  diffusion  coefficient  can  be  written  as  [80]: 


11  1 

°A(eff)  DAB(eff)  DAk(eff ) 


(14) 


Chapman  and  Enskog  is  the  most  common  method  for  theo¬ 
retical  estimation  of  gases  penetration  in  the  porous  electrodes 
while  Fuller  correlation  is  derived  based  on  the  experimental  and 
empirical  data  [81,82]. 

Table  3  shows  some  researchers  that  focused  on  mass  transfer 
in  their  model. 


3. 1 .2.  Heat  transfer  model 

In  all  the  models  that  has  been  considered  the  thermal  model, 
showed  that  the  material  properties,  chemical  kinetics,  and 
transport  properties  of  materials  used  in  a  single  cell  and  stack 
depends  on  the  temperature.  Therefore,  it  is  essential  to  develop 
a  heat  transfer  model  that  could  account  for  various  heat  effects 
in  both  of  the  solid  structure  and  the  flow  channels.  However, 
many  researchers  assumed  an  isothermal  system  in  their  model 
[38,60-62,64-66,83-87].  On  the  other  hand,  several  literatures 
studied  on  heat  transfer  modeling  in  details  in  order  to  show  the 
magnitude  of  heat  transfer  components.  Fig.  5  shows  a  comparison 
between  a  non-isothermal  and  isothermal  models.  This  figure 
depicts  the  responses  of  voltage  in  both  isothermal  and  non- 
isothermal  conditions  by  changing  a  fuel  cell  load.  As  it  can  be  seen 
in  Fig.  4,  at  the  non-isothermal  condition  the  voltage  response  is 
slow  in  comparison  with  the  isothermal  condition,  since  the  heat 
capacity  of  solid  materials  is  large  and  causing  a  slow  response. 
This  figure  shows  that  temperature  strikingly  affects  the  fuel  cell 
system  and  thus  it  is  noteworthy  to  account  energy  balance  in 
the  SOFCs  modeling.  Table  4  depicts  some  authors  that  focused  on 
heat  transfer  in  their  model. 


DAk(eff)  =  DAk 


(i)  An  energy  balance  for  the  TSOFC  inside  the  injection  tube  (SSI ) 
has  the  form: 


where  f  is  the  tortuasity. 

Using  the  Chapman-Enskog  theory  of  prediction,  the  binary 
ordinary  diffusion  coefficient  in  the  gas  phase  [73]  can  be 
calculated  as  follows  which  was  used  by  Yakabe  et  al.  [21], 
Chan  et  al.  [30],  Ni  et  al.  [41]  and  many  other  researchers 
[21 ,30,39-41 ,45,46,54,74-76]. 


DAB  =  0.0018583 


1  1  \1/2  T3/2 

Ma  +  Mb)  Po2abQdab 


(10) 


where  p  is  the  total  pressure,  crAB  =  ( oA  +  <tb)/2  is  the  collision  diam¬ 
eter  and  P2dab  is  the  collision  integral  based  on  the  Lennard-Jones 
potential  which  can  be  obtained  from  eAB  is  the  energy  of  molecular 
interaction: 


eAB 

k 


(ii) 


dp'wHmj 

r  air  air 

dt 


JW 


V  A  p'yu 

rair  air 


mj  j_jinj  ,  ninj 


air 


+  0 


onv 


(15) 


This  equation  represents  the  transport  of  energy  in  the  injec¬ 
tion  tube  by  bulk  fluid  flows  (VZ\p'nj  u'nj  H,nj  ),  convection  heat 

J  K  rair  air  air'' 

transfer  (Q^nv)  and  radiation  between  air  flow  components 

(O 

Temperature  of  the  inlet  air  stream  into  the  SSI  has  the 
strongest  effect  on  the  cell  performance  [22,43]. 

Convection  heat  transfer  in  SSI  for  TSOFC  is  between  air 
flow  temperature  and  inside  the  wall  of  injection  tube  which 
is  given  by: 


(&  =  hiti{Tit.  - 1%) 


(16) 


where  h  is  the  convection  heat  transfer  coefficient,  T1^  is  the 
air  flow  temperature  inside  the  injection  tube.  Subscript  if  is 
inside  the  of  injection  tube. 
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Table  3 

Selected  studies  concerning  on  mass  transfer  through  the  SOFCs. 


1899 


Author 

Cell  type 

Objective 

Method 

Short  coming 

Assumption 

Izzo  Jr  et  al.  [38] 

Micro-TSOFC 

-  Studying  the  gas  flow 

Gas  diffusion  was 

The  effect  of  DGM  on  fuel 

1 D,  Steady  state, 

through  an 

modeled  using  the 

cell  performance  was  not 

isothermal 

anode-supported  tubular 

dusty-gas  equations. 

compared  with  FL  and 

SOFC  and  the  subsequent 
diffusion  of  gas  through 
its  porous  anode. 

-predict  the  polarization 

Mercury  intrusion 

SMM. 

of  the  fuel  cell 

porosimetry  (MIP)  was 
used  to  experimentally 
determine 
micro-structural 
parameters. 

Bhattacharyya 

TSOFC 

-Study  the  effect  of 

In  the  process  of 

-The  effect  of  Three 

2D  diffusion, 

et  al.  [76] 

Knudson  diffusion 

validation,  phenomena 

diffusional  models  was 

isothermal  dynamic 

that  affect  the  transient 

not  considered. 

model  of  an 

response  of  the  cell 

anode-supported 

significantly  were 
identified. 

tubular  counter-flow 

-Study  the  system 

-Further  modifications 

response  and  compare 

may  be  needed  in  order 

with  experimental 

to  have  a  validated 

results. 

steady  state  and  dynamic 
model  in  a  broader 
operating  range. 

Hussain  et  al. 

Anode-supported 

-Develop  a  mathematical 

A  modified 

-It  was  not  clear  among 

1 D  steady  state, 

[37,83] 

SOFC 

model  to  describe  the 

Stefan-Maxwell 

three  diffuisional  models 

Isothermal 

transport  of 

equations  incorporating 

why  SMM  was  used  in 

multi-component  species 

Knudsen  diffusion  were 

this  model. 

inside  the  porous  anode 

used  to  model  the 

including  reaction  zone 

multi-component 

layer. 

diffusion  inside  the 
porous  anode. 

-To  anticipate  the 

-The  momentum  of 

distribution  of 

gaseous  species  in 

multi-component  species 

porous  anode  was 

in  the  electrode  and 

assumed  to  be  negligible. 

reaction  zone  layers  at 
different  current  density 
by  the  composition  of 
ethanol  reformate  fuel. 
-Study  the  effect  of  shift 

-The  effect  of  viscose  and 

reaction  and  finite 

gravity  for  ideal  gases 

reaction  zone  layer  on 
the  concentration 
overpotential. 

was  ignored. 

Lehnert  et  al.  [65] 

SOFC 

To  describe  the  transport 

Mass  transport 

-It  was  not  validated 

ID,  isothermal  steady 

of  gases  inside  the  SOFC 

phenomena  in  the  porous 

with  experimental  data. 

state 

anode  due  to  diffusion 

cermet,  such  as  diffusion 
and  permeation, 
accompanied  by 
chemical  and 
electrochemical 
reactions  internal 
reforming  and 
electrochemistry  were 
studied  and  the  resulting 
model  is  implemented  in 
a  numerical  simulation 

program. 

Suwanwarangkul 

TSOFC 

To  evaluate  the  various 

Three  types  of  models 

-  The  total  pressure 

1 D  steady  state, 

et  al.  [52] 

mass  transport  models 

including  Picks’  model, 

gradient  inside  the 

Isothermal 

by  validating  them  with 

the  dusty  gas  model  and 

electrode  was  assumed 

the  measured 

the  Stephan-Maxwell 

to  be  negligible. 

concentration 

model  inside  the  porous 

overpotential  data 

SOFC  anode  were 

obtained  from  Yakabe 

developed  to  predict  the 

et  al.  [42]. 

concentration 

overpotential. 

(ii)  An  energy  balance  for  the  PSOFC  inside  the  cathode  channel 
(SSI)  has  the  form: 


dpch'cHch'c 

•_ air  air 

dt 


vv  a  ch,c  ch,c f_rch,c  nch,c  nch,c  nch,c 

=  WAPair  Uair  Hair  +  Vconv  +  <4,d  +  <4, 


tiff 


(17) 


This  equation  represents  the  transport  of  energy  in  the  cath¬ 
ode  channel  by  bulk  fluid  flows  (V A plwulw Hinj),  convection  heat 

J  K  rair  air  avrJ' 

transfer  (Qconv),  radiation  between  air  flow  components,  cathode 
side  interconnector  wall  and  cathode  wall  (Qrad)  and  diffusional 
heat  transfer  by  oxygen  {Qdiff)- 
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Fig.  3.  Division  of  the  single  TSOFC  into  five  subsystems. 


Fuel 


SS3 


Fuel 


Anode  side  Interconnector- 


SS4 


Cathode  side  Interconnector 


SS3 


At  the  anode  surface  |  CH4+H2O 

C0+H20 


SS2 


CO+3H2  (Reforming) 
C02+  H2  (Shift) 


Fig.  4.  Division  of  the  single  DIR-PSOFC  into  five  subsystems. 


Table  4 

Selected  studies  concerning  on  heat  transfer  through  Solid  Oxide  fuel. 


Author 

Cell  type 

Objective 

Method 

Short  coming 

Assumption  software 

Inui  et  al.  [89] 

PSOFC 

Numerically  optimize  the 
operating  parameters  of 
the  cell  and  the 
temperature  of  the  inlet 
gas 

Concerning  the  calculation  of 
the  gas  flow  in  each  channel, 
the  one  dimensional  steady 
state  mass  and  energy 
conservation  equations  were 
used,  and  the  flow  was 
determined  by  integrating 
these  equations  from  the  cell 
inlet  to  its  outlet. 

The  radiation  heat  transfer 
was  not  considered 

ID  transient 

Daun  et  al.  [90] 

PSOFC 

Show  the  importance  of 
radiation  thermal  transfer. 

Two  different  radiation  solvers 
were  employed  to  solve  for  the 
radiation  source  term  in  the 
electrolyte:  the  Schuster 
Schwarzschild  two-flux 
method  [138]  and  a 
collision-based  Monte  Carlo 
method  [139]. 

-The  flow  properties  such 
as  heat  capacity,  viscosity 
were  not  considered, 

-The  enthalpy  and  entropy 
were  assumed  constant 
and  are  not  a  function  of 
temperature 

-Diffusional  heat  transfer 
was  not  considered 

2D  steady  state 

Tanaka  et  al.  [91  ] 

PSOFC 

Examine  the  effect  of  the 
heat  transfer  due  to 
radiation  on  the  cell 
operation, 

It  was  developed  a  three 
dimensional  simulation  code  of 
the  planar  SOFC  stack 

The  time  dependent 
calculation  was  performed 
only  for  the  cell 
temperature  distribution. 

3D  steady  state 

Damm  et  al.  [92] 

TSOFC 

Understanding,  predicting, 
and  quantifying  the  effects 
of  radiation  in  SOFC 
materials  and  systems. 

Modeling  surface-to-surface 
radiative  exchange,  Radiative 
transfer  in  the  flow  Channels 
and  Heat  loss  from  the  edges 

Shift/  reforming  heating 
and  diffusional  heat 
transfer  were  neglected. 

1 D  steady  state 

Brus  et  al.  [93] 

TSOFC-IIR 

To  describe  numerical 
modeling  of  the  radiative 
heat  transfer  process 

Heat  transfer  between  the  fuel 
reformer  surface  and  all  other 
surfaces  facing  the  reformer 
surfaces  was  modeled 

-The  flow  properties  such 
as  heat  capacity,  viscosity 
were  not  considered 

-The  enthalpies  were 
assumed  to  be  constant 

STanchez  et  al.  [58] 

TSOFC 

To  simulate  a  complicated 
heat  transfer  model  and 
compare  the  errors  used 
for  those  with  simple 
models  with  the  complex 
model. 

The  work  was  divided  in  two 
parts.  First,  two  models  for 
describing  convective  heat 
transfer  was  proposed.  For  the 
radiation  heat  transfer,  two 
models  were  presented  again 

-The  velocity  and  mass 
flow  were  assumed 
constant  along  the  length 
therefore  the  conduction 
heat  transfer  coefficient 

was  constant. 

1 D  steady  state 
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Convection  heat  transfer  in  SSI  for  PSOFC  is  between  air  flow 
and  both  cathode  channel  and  cathode  side  interconnector  walls 
which  is  given  by: 

Q^onv  =  hint^T-nli  ~  Tair )  +  W  Wo  “  T air )  (18) 

where  subscripts  int,q  and  ct,c0  are  inside  the  cathode  side  inter¬ 
connector  and  outside  the  cathode  side  of  the  cell  tube,  respectively, 
and  superscript  cat  is  the  cathode  side. 

An  empirical  Nusselt  number  correlation  is  used  to  calculate  the 
convection  heat  transfer  coefficient: 


NU  =  f(Re,  Pr)  (19) 

NU  =  ^  (20) 

D  is  the  fuel  cell  diameter. 

By  Hayes  and  Kolaczkowski  [88]  expression  Nusselt  number 
yields  as  the  follow: 


NU  =  3.095  +  8.933 


lOOOV0'5386 

-qT)  exp 


6.7275  \ 
Gz  ) 


(21) 


where  the  Greatz  number  Gz  is  given  by 

Gz  =  -RePr  (22) 

X 

x  is  the  axial  position,  Re  and  Pr  are  the  Reynolds  number  and  Prantl 
number,  respectively. 

For  a  fully  developed  laminar  flow  approximation  at  constant 
wall  temperature  Nusslet  number  is  assumed  to  be  3.66. 

For  Re  >  2300,  0.5  <Pr  <  2000  and  L>  Dh  Gienelinski’s  equation  is 
applicable: 


NU  = 


Pr(Re  -  1000)(//8) 
l  +  12.7(Pr2/3-l)x/(78 


(23) 


/  =  ( 0.79  ln(Re)  -  1.64)-1 


(24) 


Table  5 

Different  values  of  conduction/convection  heat  transfer  coefficient. 


Authors 

hc  (kj/m2  sK) 

I<  (kj/m  s  K) 

Aguiar,  P.  et  al.  [95] 

hC}a  =  hc,c  =  0.1 

2  x  10-3 

Damm  et  al.  [92] 

bc,a  =  bc,c  =  0.6 

I<a  =0.00584 

Kc  =  0.001 86 

Keiec  =  0.00216 

Kair  channel  =  0.000047 

Kfuel  channel  =  0.0002 

Denver  et  al.  [92] 

hc,a  =  hCfC  =  0.302 

Xue  et  al.  [84] 

hc,a  =  2.987 
hc,c  =  1.322 

I<a=Kc  =  Kelec  =  0.00253 

For  PSOFC  in  SSI,  oxygen  diffuses  into  the  cathode  in  order  to 
participate  in  electrochemical  reaction,  thus  the  diffusional  heat 
transfer  leads  to: 

Qdl ff  =  VNo2Ho2  (27) 

Air  mainly  consists  of  non-polar  nitrogen  and  oxygen  molecules 
and  therefore,  is  non-interaction  (transparent)  with  radiation  heat 
transfer  at  the  moderate  temperatures  and  pressures  within  SOFCs 
[93].  Therefore,  no  radiation  heat  transfer  occurs  between  the  air 
properties.  So  in  SSI  for  both  TSOFC  and  PSOFC  no  radiation  heat 
transfer  occurs. 

3.1.3.  Momentum  model 

Velocity  affects  the  film  heat  transfer  coefficient  on  the  anode 
and  cathode  side.  An  increase  in  the  gas  flow  velocity  leads  to  an 
increase  in  the  film  heat-transfer  coefficient,  causing  the  removal 
of  more  energy  from  the  cell  tube,  resulting  in  a  decrease  in  the 
cell-tube  temperature  [22].  On  the  other  hand,  velocity  is  affected 
by  reactions,  mass  transfer,  density  change,  etc.  The  flow  velocities 
can  be  modeled  by  the  momentum  conservation. 


where  Dh  is  hydraulic  diameter,  L  the  length  of  fuel  cell  and/friction 
factor.  However  this  equation  has  complexity  for  calculating  the 
convection  coefficient. 

Colburn’s  equations  is  easier  to  use  but  have  errors  around  20% 
[58]: 


(i)  A  momentum  balance  inside  the  TSOFC  inside  the  injection  tube 
(SSI)  based  on  Navier-Stokes  equation  leads  to  [94]: 


VPinj  +  Vrinj  +  Finj 


air 


air 


air 


(28) 


For  laminar  flow: 

NU  =  0.332 Pr1/3Re1/2 

(25) 

For  turbulent  flow: 

NU  =  0.023Pr1/3Re4/5 

(26) 

In  some  studies  constant  values  for  h  and  I<  at  transient  condi¬ 
tion  are  assumed  [83,89].  In  the  study  of  D.  S’anchez  et  al.  [58]  the 
effect  of  constant  and  variable  convective  heat  transfer  on  fuel  cell 
was  compared.  The  results  showed  that  these  two  models  substan¬ 
tially  affect  the  temperature  at  both  ends  of  the  fuel  cell  tube  while 
there  is  no  impact  on  the  shape  of  temperature  curve.  Some  authors 
like  Hajimolana  and  Soroush  [22],  Ji  et  al.  [29]  and  Qi  et  al.  [43]  cal¬ 
culated  thermal  conductivity,  convection  heat  transfer  coefficient, 
viscosity,  specific  heat  capacity,  enthalpy  and  specific  resistivity  as 
a  function  of  temperature.  The  values  of  conduction  and  convection 
heat  transfer  coefficient  existing  in  the  literatures  are  illustrated  in 
Table  5. 

Energy  transportation  inside  and  outside  of  the  mass  media 
occurs  by  diffusion  of  components.  Diffusional  heat  transfer  is 
available  in  [22,39,41,43,46,48,90-92].  However,  this  type  of  heat 
transfer  was  not  involved  in  some  literatures  because  it  was 
assumed  that  it  is  not  magnitude  compared  to  the  conduction  heat 
transfer  in  porous  media  [61,67]. 


where  (  V  puu)  is  the  momentum  rate  flowing  in  and  out  of  the 
SSI,  VP  the  pressure  gradient,  r  is  the  shear  stress,  and  F  is 
prescribed  body  force  (i.e.  gravity  force). 

(ii)  A  momentum  balance  for  the  PSOFC  inside  the  cathode  channel 
(SSI)  based  on  Navier-Stokes  equation  leads  to  [94]: 


9(  pcKc  uch’c) 

Krair  air  J 


dt 


+  {V  pch’cuch’cuch,c) 

K  rair  air  air  J 


_L  yrch,c  pch,c 
r  air  '  air 


—N  o2  Mq2u 


ch,c 

air 


ypch,c 

air 


(29) 


where  N02Mq2u  is  the  momentum  out  flow  rate  resulting  from 
02  participation  in  the  electrochemical  reaction. 


Ferguson  et  al.  [50],  Janardhanan  et  al.  [67]  and  other 
researchers  [50,67,78,95-97]  neglected  the  momentum  con¬ 
servation  in  their  model,  while  many  researchers  like  Haji¬ 
molana  and  Soroush  [22],  Bove  et  al.  [32],  Autissier  et  al. 
[35],  Qi  et  al.  [43],  Suwanwarangkul  et  al.  [51]  and  others 
[32,35,39,43,45,51,53,54,61,63,66,69,77,78,83,86,91,98-103]  con¬ 
sidered  the  momentum  conservation  in  their  model.  Bhattacharyya 
et  al.  [54]  compared  the  models  with  and  without  momentum  con¬ 
servation  and  found  that  the  change  in  velocity,  especially  in  the 
cathode  side  channel  is  significant.  When  momentum  conserva¬ 
tion  is  involved  in  the  model  the  concentration  of  oxygen  always 
remains  higher,  hence  the  oxygen  concentration  is  also  higher  at 
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the  TPB.  It  was  observed  that  momentum  conservation  is  important 
for  high  power  range. 

The  pressure  gradient  VP  can  be  approximately  determined 
using  Poiseuille’s  Law  [104]. 

By  calculating  the  friction  factor/the  wall  shear  stress  yields: 


3.2.  Model  equations  for  subsystem  2  (SS2 ) 

This  subsystem  for  TSOFC  includes  the  solid  of  injection  tube 
and  for  PSOFC  is  the  solid  of  the  cathode  side  interconnector.  In 
this  subsystem  mass  and  energy  balances  can  be  accounted. 


/ 


2r 

pu2 


(30) 


/is  constant  for  laminar  flow  and  for  turbulent  flow  it  depends  on 
Reynolds  number  and  surface  roughness  [81]. 

In  the  case  of  the  channel  volume  that  is  represented  by  porous 
medium  (for  the  case  of  effective  parameters)  Darcy  equation  can 
describe  the  balance  between  the  force  from  the  pressure  gradient 
VP  and  the  frictional  resistance  s  from  the  solid  material  which  is 
given  by: 


-VP  = 


irUp 


(31) 


where  Up  is  the  pore  velocity  vector,  s  is  the  porosity  (the  volume 
fraction  of  void  space),  pg  is  the  dynamic  viscosity  of  the  gas  mix¬ 
ture  and  K  is  the  permeability  which  depends  on  micro-structure 
of  the  porous  electrodes. 

For  random  packing  systems  of  spherical  particles,  Caman- 
Kozeny  correlation  is  generally  used  to  evaluate  the  flow 
permeability  as: 

£3 

I<  = - - — —  (32) 

kK(  1  -£)  A2 

where  the  Kozeny  constant  kk  is  about  5  for  porous  media  made  of 
spherical  particles. 

Aq  is  the  specific  area  based  on  the  solid  volume: 


3.2. 1 .  Heat  transfer  model 

(i)  An  energy  balance  for  TSOFC  in  the  injection  tube  (SS2)  leads 
to: 

dC  T 

=  (36) 

where  m  is  the  mass  and  Cp  is  the  heat  capacity. 

Convection  heat  transfer  in  SS2  for  TSOFC  is  between  inside 
injection  tube  and  air  flow  and  also  between  outside  the  injec¬ 
tion  tube  and  the  air  flow  inside  the  SS3  and  can  be  written  as 
follows: 

Qllnv  =  MC  -  V  +  MC  -  Tito)  (37) 

Heat  transfers  inside  the  gases  and  through  the  elec¬ 
trodes,  electrolyte  and  interconnectors  walls  with  conduction 
heat  transfer,  however,  the  thermal  conductivity  of  gases  is 
mainly  very  low  (0.025  W/mK)  [66].  Conduction  heat  transfer 
described  by  Fourier’s  equation  which  is  given  by: 

Qcond  =  VKVT  (38) 

I<  is  the  thermal  conductivity  of  the  solid  structure. 

Subscript  it0  is  the  outside  of  the  injection  tube. 

Thermal  radiation  within  tubular  SOFC  in  SS2  occurs  between 
injection  tube  and  inside  the  cell  tube  walls  which  is  given  by 
eq  (38): 


6  nei  +  ( 1  -  nei)a2 
del  ne/  +  (l  -  nei)a3 


(33) 


a  is  bi-model  size  distribution. 

The  volume  flow  rate  per  unit  area  through  the  porous  structure 
is  obtained  by  multiplying  Up.  The  pore  Reynolds  number  is  defined 
by: 


=  08') 
lvrad 

Rrad  is  the  radiation  heat  transfer  resistance  which  is  given  by 
[80]: 


1  ~  £ct  1  1  ~  £it  At0 

Set  Fct-it  £it  Act{ 


(39) 


Pg 


U, 


p 


lp 


pg 


(34) 


where  pg  is  the  gas  mixture  density.  Eq.  (34)  is  valid  for  Rep  <  1  [105] 
when  pressure  loss  is  dominated  by  pore  friction  and  the  inertia 
terms  usually  present  in  a  momentum  equation  are  insignificant. 
The  Darcy  equation  describes  the  flow  in  the  porous  structure  well 
away  from  the  walls.  It  cannot  model  a  no-slip  condition  at  a  wall 
nor  the  resulting  boundary  layers.  It  is  found  that,  although  the 
boundary  layer  thicknesses  are  typically  only  about  1-3%  of  the 
ceramic  support  thickness,  Eq.  (31)  can  be  modified  to  allow  their 
calculation.  This  is  done  by  adding  the  so-called  ‘Brinkman’  term  to 
the  Darcy  equation  which  then  becomes: 


-VP  = 


^~§—Up  —  V.(/i,gV)[/p 


(35) 


where  the  normal  stresses  in  the  final  term  as  usual  are 
ignored.  Darcy  equation  was  used  by  several  researchers  like 
Yakabe  et  al.  [21],  Bove  et  al.  [32],  Ho  et  al.  [46]  and  others 
[21,32,38,39,46,64,69,79,101-103,106,107]. 

Pressure  drop  in  both  anode  and  cathode  sides  was  neglected 
by  some  authors  like  Shi  et  al.  [85]  Aguiar  et  al.  [108,109]  and 
others  [78,85,96,97,108,1 10-1 13].  However,  the  temperature  vari¬ 
ation  along  the  cell  can  be  accurately  predicted  by  considering  the 
pressure  gradient  among  the  fuel  cell  [114].  At  the  cathode  side, 
neglecting  the  pressure  gradient  can  lead  to  make  an  error  in  cal¬ 
culating  of  concentration  overpotential  by  about  10%  under  typical 
working  conditions  [66]. 


(ii)  An  energy  balance  for  PSOFC  in  the  cathode  side  interconnector 
(SS2) leads  to: 


m 


dCPint,Ji 


int,c 


int,c 


dt 


f)int,c  nint,c  nint,c 
<Zcony  "I"  40nd  ^  4 ad 


(40) 


Convection  heat  transfer  for  PSOFC  is  between  the  solid  of  the 
air  side  interconnector  and  air  flow  in  SS2  which  is  given  by  Eq. 
(41): 


Qeorw  =  hi„t,Ci(Tai 


air 


(41) 


Thermal  radiation  within  PSOFC  in  SS2  occurs  between  air 
side  interconnector  and  cathode  wall  as  follows: 


O  int,c  &  / y 4  _ rrcaH  \ 

4-ad  —  P  ,  ^  ct,c0  1  int,  z  ' 
1Krad 


(42) 


1  —  Scat  1 

Scat  Feat- int 


+ 


1 


£jnt 

4nt 


Ant 

Acat 


(43) 


Because  of  the  high  operating  temperatures  of  SOFCs  (typically 
800-1200 1<),  there  has  been  some  estimation  that  thermal  radi¬ 
ation  may  be  an  important  factor  of  heat  transfer  within  them. 
Radiation  heat  transfer  has  been  involved  by  Hajimolana  and 
Soroush  [22],  Damm  et  al.  [89],  Aguiar  et  al.  [108],  Iora  et  al. 
[91],  Jia  et  al.  [115],  Achenbach  [116]  and  other  several  mod¬ 
els  in  [26-28,35,44,51,61,75,77,83,89,91,96,97,108,110,115-120]. 
Detailed  discussions  of  radiation  heat  transfer  within  SOFCs  are 
given  by  Hajimolana  and  Soroush  [22],  Qi  et  al.  [43],  Ho  et  al.  [46], 
Sanchez  et  al.  [58],  Damm  et  al.  [89],  Brus  et  al.  [121],  Wang  et  al. 
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[98],  Calise  et  al.  [122],  Haynes  [123].  However,  some  researchers 
ignored  radiation  heat  transfer  due  to  considering  a  uniform  tem¬ 
perature  everywhere  in  the  cell  [55],  or  because  of  assuming  a  thin 
wall  of  channels  [124].  In  Iwai,  H.,  etal.  [125]  model  it  was  not  clear 
why  the  radiation  heat  transfer  was  not  tacked  into  account.  Wang 
et  al.  [126]  and  Yakabe  et  al.  [31]  neglected  the  effect  of  radiation 
heat  transfer  in  their  model  due  to  a  very  small  amount  compared 
to  other  types  of  heat  transfer. 

In  fact,  the  magnitude  of  thermal  radiation  depends  on  two 
parameters,  temperature  and  the  thickness  of  electrodes  and  elec¬ 
trolyte.  At  a  greater  temperature  and  a  thicker  of  the  layers,  the 
radiative  transfer  becomes  more  significant  as  compared  to  other 
heat  transfer  components. 

Calise  et  al.  [122]  presented  a  model  of  a  tubular  SOFC  in  which 
the  radiation  heat  transfer  was  modeled  in  details.  It  was  observed 
that  radiation  heat  transfer  has  a  strong  effect  on  temperature  dis¬ 
tribution  among  the  fuel  cell  (about  70%).  The  model  was  developed 
by  removing  some  simplification  used  in  the  study  of  Companari 
[127]  and  Stiller  groups  [28].  In  Companari  [127]  the  radiation 
heat  transfer  between  air  injection  tube  and  the  SOFC  tube  was 
neglected  while  F.  Calise  group  [122]  considered  that  the  radia¬ 
tive  heat  transfer  is  predominately  through  the  radial  distribution 
among  the  fuel  cell  tube  and  air  injection  and  is  strikingly  impor¬ 
tant. 

Sanchez  et  al.  [58]  performed  a  steady-state  model  of  a  tubular 
SOFC  including  convection  and  radiation  heat  transfers.  Two  mod¬ 
els  for  convection  heat  transfer  were  studied.  In  the  first  model 
film  coefficient  was  a  function  of  local  temperature  and  compo¬ 
sition.  The  local  value  for  the  heat  transfer  coefficient  and  the 
thermal  entry  length  of  the  cell  is  clear  with  this  model.  In  the 
second  model,  an  average  value  of  the  heat  transfer  coefficient 
for  the  whole  length  of  the  fuel  cell  channel  was  assumed.  Based 
on  this  study,  it  was  observed  that  a  simple  model  can  be  used 
to  evaluate  the  fuel  cell  performance  with  satisfactory  accuracy, 
although  there  is  a  need  to  calculate  local  temperatures.  For  radi¬ 
ation  heat  transfer,  two  models  were  performed  again.  The  first 
model  accounted  the  radial  radiation  exclusively  and  thus,  radiative 
exchange  between  adjacent  cells  is  neglected.  The  second  model 
considered  the  radiation  completely  in  all  directions.  This  causes 
to  increase  significantly  the  complexity  of  model.  It  was  found  that 
a  simple  model  based  on  a  constant  film  convective  heat  transfer 
coefficient  and  a  simple  radiation  model  is  enough  for  predicting 
the  fuel  cell  performance  without  internal  information  requiring 
and  using  a  complicate  model  is  not  necessary.  A  complete  heat 
transfer  model,  however,  is  necessary  for  evaluating  the  internal 
fuel  cell  temperature  accurately. 

Conduction  heat  transfer  was  tacked  into  account  by 
Stiller  et  al.  [28],  Negata  et  al.  [44],  Ho  et  al.  [45,46],  Fer¬ 
guson  et  al.  [50],  Sanchez  et  al.  [58]  and  other  researchers 
[28,44-46,50,58,63,67,75,77,78,89,92,99,101,102,108,110,115,117, 
128,129]  in  order  to  model  the  distribution  temperature  of  fuel 
cell  more  accurately.  However;  no  conduction  heat  transfer  was 
used  by  Hajimolana  and  Soroush  [22],  Costamagna  et  al.  [26], 
Burt  et  al.  [27],  Qi  et  al.  [43]  and  Achenbach  [116].  In  this  studies 
conduction  heat  transfer  was  neglected  due  to  assuming  a  thin 
thickness  of  electrodes,  electrolyte  and  interconnectors  walls  or  to 
make  simple  the  heat  transfer  modeling  for  calculating.  It  was  also 
found  that  at  low  temperature  the  conductivity  through  the  YSZ 
electrolyte  becomes  limiting  [76]. 


3.3.  Model  equations  for  subsystem  3  (SS3 ) 

This  subsystem  for  TSOFC  includes  the  air  inside  the  space 
between  the  cell  and  injection  tubes  and  for  PSOFC  includes  the 
fuel  side  interconnector. 


3.3.  \ .  Mass  transfer  model 

(i)  A  mass  balance  on  the  air  inside  SS3  for  TSOFC  takes  the  form: 


(44) 


A  mass  balance  on  oxygen  leads  to: 


+  VKi>“rty 02)  =  -VNo2 


(45) 


The  effect  of  02  diffusion  on  cell  voltage  and  current  is  not  sig¬ 
nificant  because  the  concentration  of  02  in  the  air  is  21%  [25,34]. 


3.3.2.  Heat  transfer  model 

(i)  An  energy  balance  for  TSOFC  inside  the  SS3  leads  to: 


=  VAp 


cat 

air  uair 


catj_fcat  ,  neat 


air 


+  0 
'  v<ci 


onv 


(46) 


which  accounts  for  energy  flowing  into  and  out  of  SS3  with 
air  ( V  ApcatucatHcat),  convectional  heat  transfer  and  diffusional 

v  rair  air  air 

heat  transfer  via  oxygen  into  the  cathode  (Q^S), 

Convective  heat  transfers  from  the  cell  and  injection  tubes 
are  as  follows: 


QcCoL  =  hito(Tito  -  T%)  +  hct^Tct,  -  Cf) 

(ii)  An  energy  balance  for  PSOFC  inside  the  SS3  leads  to: 


(47) 
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(48) 


Convection  heat  transfer  for  PSOFC  is  between  the  solid  of  the 
fuel  side  interconnector  and  fuel  flow  in  SS7  and  leads  to: 


int,a  ( rrano  ycfr,a\ 

1 ^conv  —  ,tint,aAiint,i  lfuel} 


(49) 


where  superscript  ano  denotes  the  anode  side. 

Thermal  radiation  within  PSOFC  in  SS3  occurs  between  fuel 
side  interconnector  and  anode  side  walls  and  is  as  follows: 


r\  int,a  _  &  (t4  _ rrano^  \ 

^rad  —  P  .  ^  ct,a0  Jint,i  ' 
lvrad 


(50) 


where  subscript  ct,a0  denotes  the  outside  the  anode  side  of  the 
cell  tube. 


3.3.3.  Momentum  model 

Momentum  conservation  on  the  air  inside  SS3  for  TSOFC  takes 
the  form: 


T  (V  Pair^ 
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(51) 


3.4.  Model  equations  for  subsystem  4  (SS4) 

This  subsystem  for  both  TSOFC  and  PSOFC  is  the  cell  tube  includ¬ 
ing  electrodes  and  electrolyte.  Electrochemical  reaction,  reforming 
and  shift  reactions  occur  in  this  subsystem. 

3.4 A.  Electrochemical  reaction 

High  SOFC  performance  depends  on  optimal  electrochemical 
reactions  and  mass  transport  processes.  Since  experimental  stud¬ 
ies  on  SOFC  are  expensive,  time-consuming  and  labor-intensive, 
quantitative  mechanistic  models  for  the  cell  PEN  structure  are 
essential  for  SOFC  technology  development.  Electrochemical  reac¬ 
tion  occurs  inside  the  fuel  cell  which  as  a  result  produces  voltage 
and  current.  At  cathode  side  oxygen  ions  (with  a  negative  charge) 
migrate  through  the  crystal  lattice.  The  oxygen  is  supplied,  usu¬ 
ally  from  air,  at  the  cathode.  Hydrogen  and/or  carbon  monoxide 
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(CO)  in  the  fuel  stream  reacts  with  oxide  ions  (02~)  from  the 
electrolyte  to  produce  water  or  C02  and  deposit  electrons  into 
the  anode.  The  electrons  pass  outside  the  fuel  cell,  through  the 
load,  and  back  to  the  cathode,  where  oxygen  from  the  air  receives 
the  electrons  and  is  converted  to  oxide  ions,  which  are  then 
injected  into  the  electrolyte.  Generating  efficiencies  can  range  up 
to  ~60%.  In  one  configuration,  the  SOFC  consists  of  an  array  of 
tubes.  Another  variation  includes  a  more  conventional  stack  of 
disks.  The  high  operating  temperature  of  the  SOFC  offers  the  pos¬ 
sibility  of  internal  reforming.  As  in  the  molten  carbon  fuel  cells, 
CO  does  not  act  as  a  poison  and  can  be  used  directly  as  a  fuel. 

{H2  +02“  -*  H20  +  2e 
CO  +  O2-  ->  C02  +  2e 
CH4  +  402-  -*  2H20  +  C02  +  8e 

Cathode  Side  Reaction  {02(g)  +  4e  ->  202- 

The  1000°C  operating  temperature  of  SOFCs  requires  a  signif¬ 
icant  start-up  time.  The  cell  performance  is  very  sensitive  to  the 
operating  temperature.  A  10%  decrease  in  temperature  can  lead  to 
a  12%  drop  in  cell  performance,  because  of  a  higher  internal  resis¬ 
tance  to  the  flow  of  oxygen  ions  [46].  The  high  temperature  also 
demands  that  the  system  include  significant  thermal  shielding  to 
protect  personnel  and  to  retain  heat.  Such  requirements  are  accept¬ 
able  in  utility  applications  but  not  in  small  portable,  transportable, 
or  most  transportation  applications. 

It  is  significant  that  a  SOFC  can  use  CO,  as  well  as  hydrogen 
as  its  direct  fuel.  When  both  species  of  H2  and  CO  are  present  in 
the  system,  CO  mostly  participates  in  the  water-gas  shift  reac¬ 
tion  of  equation  [94]  rather  than  in  the  electrochemical  process 
[130].  Flowever,  current  produces  by  the  sum  of  hydrogen  and  car¬ 
bon  monoxide  oxidation  at  the  anode  side  by  the  electrochemica 
reaction.  The  rate  of  H2  oxidation  is  2-3  times  more  than  that  of 
CO  oxidation  depending  on  the  operationg  temperature  [131].  Ho 
et  al.  [132]  assumed  that  the  rate  of  CO  oxidation  is  3  times  less 
that  of  H2  oxidation  due  to  high  operating  temperature  of  the  fuel 
cell.  On  the  other  hand,  Hofmann  et  al.  [130]  assumed  that  only 
H2  participate  at  electrochemical  reaction  due  to  its  fast  kinetic 
reaction. 

It  is  noteworthy  that  when  the  shift  reaction  is  at  equilib¬ 
rium  condition,  CO  and  H2  oxidation  produces  the  same  Nernst 
potential  in  the  fuel  cell  [131].  Many  authors  like  Cayan  et  al. 
[40],  Yakabe  et  al.  [42],  Qi  et  al.  [43]  and  other  researchers  in 
[40,42,43,61,68,69,112,116,133-135]  considered  CO  oxidation  in 
their  model  to  make  their  model  more  realistic.  However,  some 
researchers  like  Costamagna  et  al.  [26],  Stiller  et  al.  [28],  Ho  et  al. 
[46]  and  also  [26,28,46,102,103,106,110,113,117,120,136-138] 
ignored  occurrence  of  CO  oxidation  at  the  anode  side  due  to  domi¬ 
nating  of  H2  over  CO  in  the  charge  transfer  chemistry  It  is  estimated 
that  about  98%  of  current  is  produced  by  H2  oxidation  in  common 
situations  and  therefore  CO  seems  to  play  a  minor  role  and  can 
safely  be  neglected  [139]. 

Aloui  and  Halouani  [69]  compared  the  fuel  cell  performance 
with  and  without  carbon  monoxide  oxidation  at  the  anode  side.  It 
was  observed  by  L.  Petruzzi  et  al.  [140]  that  CO  oxidation  at  the  elec¬ 
trochemical  reaction  influences  the  cell  efficiency,  fuel  utilization, 
thermodynamic  efficiency  and  electric  power.  Actually,  oxidation 
temperature  heavily  depends  on  CO  and  as  a  result,  CO  oxidation 
can  affect  the  cell  performance. 

Methane  can  also  possibly  be  electrochemically  oxidized  in  the 
case  of  absence  or  near-absence  of  steam  feed  [93,104].  How¬ 
ever,  with  high  content  of  steam  the  fast  reforming  reaction 
will  dominate  over  any  oxidation  of  methane  and  can  safely  be 
neglected. 


Incidentally,  a  pure  voltage  source  is  usually  referred  to  as 
an  electromotive  force  (emf).  The  electromotive  force,  reversible 
open-circuit  cell  voltage  (denoted  by  U0cv )  is  given  by  the  Nernst 
equation: 


RT 

Uqcv  =  U°(T)  +  —  In 


xH,,ax°;5(Pc/Pre/) 


0.5 


*H20,a 


(52) 


Equation  [52]  is  hold  for  hydrogen,  where  U0cv  is  the  standard 
potential  cell,  R  is  the  universal  gas  constant,  F  is  Farady  constant, 
P  is  the  pressure,  H°(T)  is  the  standard  cell  potential,  and  x  is  the 
mole  fraction.  If  the  shift  reaction  is  assumed  to  be  at  equilibrium 
the  Uocv  of  Co  is  equal  to  that  of  CO. 

In  practice,  cell  emf  depends  on  temperature  and  concentra¬ 
tion  of  reactants  and  products.  If  the  concentration  of  reactants 
increases  relative  to  the  products,  the  cell  reaction  becomes  more 
spontaneous  and  the  emf  increases.  As  the  cell  operates,  the  reac¬ 
tants  are  used  up  as  more  product  is  formed  causing  the  emf  to 
decrease. 

For  a  fuel  including  H2  and  CO  the  ideal  potential  for  reaction 
using  the  Nernst  equation  is  as  follows: 

u»“-u°<r>+5{"H>in(^)+"coin(^) 

-»H1l0(^)-nco|n(^+"l°(^)}  (53) 

ne  is  the  total  number  of  electrons  transferred  and  given  by  [49]: 
ne  =  nH2  x  2  +  nCo  x  2  (54) 

However,  the  actual  cell  voltage  (H)  is  less  than  its  theoretical 
open  circuit  voltage  because  it  strongly  affected  by  several  irre¬ 
versible  losses  including  activation  losses  due  to  irreversibility  of 
electrochemical  reactions  at  the  three-phase  boundary,  concen¬ 
tration  losses  due  to  mass  transport  resistance  in  the  electrode, 
especially  for  thick  anodes  as  in  an  anode-supported  SOFC  [50]  and 
ohmic  losses  due  to  ionic  and  electronic  charge  transfer  resistances. 
Actual  voltage  is  given  by: 


U  —  U(DCV  ~  hohm  ~  hact  ~  hconc  (55) 

The  rates  of  consumption  of  hydrogen  and  oxygen  by  the  elec¬ 
trochemical  reactions  to  generate  an  electric  current  of  /  are  given 
by: 


The  consumption  of  the  reactants  is  accompanied  by  the  pro¬ 
duction  of  water  at  the  following  rate: 


%  =  (^p)j-  (57) 

3.4.2.  Electrochemical  losses 

3.4.2. 1.  Activation  Losses.  The  activation  polarizations  are  the 
results  of  the  kinetics  involved  with  the  electrochemical  reactions. 
It  becomes  an  important  loss  when  the  current  is  low  because  at 
low  current  density  (j)  the  reactants  must  overcome  an  energy 
barrier  named  activation  energy  (Fsct)  to  drive  the  electrochemi¬ 
cal  reactions  at  the  electrodes-electrolyte  interface  and  this  barrier 
leads  to  the  polarization.  The  activation  barrier  is  the  result  of  many 
complex  electrochemical  reaction  steps  where,  typically,  the  rate- 
limiting  step  is  responsible  for  the  polarization.  One  can  account 
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for  the  anode  and  cathode  activation  polarizations  using  the  well 
known  Butler-Volmer  equation  [30]: 


as  a  Taylor  series.  Therefore,  by  ignoring  the  terms  with  order 
higher  than  unity  the  following  equation  yields: 


j  =j  o  exp 


^  ®E)  act, exp  f  ^  ^)^lact,i^i^ 


RT 


RT 


(58) 


The  equation  explains  the  net  anodic  and  cathodic  current 
density  due  to  an  electrochemical  reaction.  Where  subscript  i  rep¬ 
resents  the  species  H2,  CO  and  02,  j  is  the  current  density,  n  is  the 
number  of  electrons  transferred  in  the  electrochemical  reactions, 
Tjact  is  the  activation  polarization,  a  is  the  transfer  coefficient  and 
its  value  is  usually  0.5,  and  j0  is  the  exchange  current.  When  f3  =  0.5, 
Eq.  (58)  can  be  expressed  as  follows: 

j  =  2j„sinh(^)  (59) 


The  total  current  density  is  given  by: 


ltot  j°2  .  (60) 

Jtot  =J  h2  +Jc  0 

At  the  equilibrium  potential  the  electrode/electrolyte  interface 
still  has  electron  transfer  processes  going  on  in  both  directions. 
The  ongoing  current  in  both  directions  at  the  dynamic  equilibrium 
electric  potential  differences  is  called  exchange  current  density. 
The  exchange  current  density  is  a  factor  to  measure  the  electro- 
catalytic  activity  of  TPB  for  a  given  electrochemical  reaction.  It  is 
also  a  significant  factor  for  calculating  the  activation  polarization. 
Exchange  current  density  depends  on  many  parameters  such  as 
concentration  of  reactants  and  products  at  the  electrodes,  temper¬ 
ature,  pressure,  the  microstructure  and  electrocatalytic  activity  of 
the  electrode,  and  even  conductivity  of  electrolyte  [104].  Exchange 
current  density  with  respect  to  reactant  concentrations  depen¬ 
dency  can  be  calculated  as  follows  [30]: 

Jo=j%nxp  (6i) 

I<=  1 


where  j'q  is  a  constant,  and  yk  and  XK  are  the  reaction  order  and  mole 
fraction  (at  the  electrode-electrolyte  interface)  for  the  Kth  species, 
respectively  [30,62].  A  high  exchange  current  density  causes  the 
reaction  to  occur  rapidly,  so  a  good  fuel  cell  performance  can  be 
expected. 

Solving  the  Butler-Volmer  equation  with  a  transfer  coefficient 
of  0.5  in  the  high-current-density  regime  (where  the  cell  typically 
operates)  leads  to: 


El  act,  a 
EJ  act,  c 


(62) 


where  jo, a  and  jo.c  are  the  anode  and  cathode  exchange  currents. 

Eact,ano  and  Eact,cat  are  energy  activation  of  anode  and  cathode, 
PTPB  is  the  pressure  of  components  inside  the  triple  phase  boundary, 
and  A  is  the  surface  area  of  the  cell. 

High  activation  polarization  and  low  activation  polarization  are 
two  possible  cases  of  Eq.  (58).  Under  high  activation  polarization, 
the  first  term  of  Eq.  (58)  will  be  much  greater  than  the  second  term, 
so  the  simplified  equation  is  given  as  follow: 


— (^Wi) 

which  is  called  the  Tafel  equation. 

At  low  activation  polarization,  the  term  ( ariactiFIRT )  in  Eq.  (58) 
will  be  much  less  than  and  the  exponential  terms  can  be  expanded 


which  is  a  relation  that  shows  linear  current  potential.  The  issue 
of  applying  the  range  of  activation  polarization  for  using  the  Tafel 
equation  or  linear  current  relation  is  discussed  by  S.H.  Chan  et  al. 
[30]. 

Activation  polarization  in  high  temperature  for  the  case  of  SOFCs 
is  usually  small  [141].  It  also  depends  on  CO  and  H2  concentra¬ 
tion  and  decreases  when  the  amount  of  these  components  in  fuel 
increases  due  to  the  amelioration  of  the  anode  exchange  current 
density. 


3.422.  Concentration  losses.  Concentration  losses  are  losses  asso¬ 
ciated  with  concentration  variation  of  the  critical  species  due  to 
mass  transport  processes.  There  are  usually  two  sources  of  loss 
that  are  due  to  mass  transport  loss:  (i)  diffusion  between  the  bulk 
flows  and  cell  surfaces,  and  (ii)  transport  of  reactants  and  prod¬ 
ucts  through  electrodes.  Therefore,  the  concentration  polarization 
is  highly  dependent  on  the  gases  used,  as  well  as  the  distance 
through  which  the  gases  must  diffuse.  Pore  volume  percentage, 
as  well  as  diffusion  length,  can  be  varied  to  optimize  these  prop¬ 
erties.  For  similar  geometries,  cathode  concentrations  are  much 
larger  than  anode  concentrations,  because  of  the  lower  diffusivi- 
ties  of  02/N2  in  the  cathode  than  H2/H20  in  the  anode.  The  anodic 
and  cathodic  concentration  losses  can  be  calculated,  respectively, 
as  the  following: 


El  cone,  a 


El  cone,  c 


RT 


1=1 


(65) 


where  X*  is  the  component  symbol  for  the  ith  species  which  may  be 
involve  in  the  fuel  or  oxidizer.  i/j  and  v'a  are  the  stoichiometric  coef¬ 
ficients  for  the  fuel  and  oxidizer  mixtures,  respectively,  and  v'' .  and 
v"  .  are  the  stoichiometric  coefficients  of  the  ith  product  species 
in  the  fuel  and  air  channels,  respectively.  [Xz  ] *  and  [Xz  ]s  denote 
the  molar  species  concentration  in  the  channels  and  triple  phase 
boundary,  respectively. 

This  type  of  polarization  loss  is  considerable  for  high  fuel  uti¬ 
lization  [142]  also  normally  reduced  at  high  current  density.  A 
dramatic  change  happens  for  an  anode-supported  tubular  SOFC 
which  requires  high  fuel  utilization  due  to  a  small  error  in  the 
concentration  overpotential  calculation. 


3.42.3.  Ohmic  losses.  Ohmic  losses  happen  due  to  existence  of 
resistance  for  ions  transferring  through  the  electrolyte  and  elec¬ 
trons  through  the  electrodes  and  current  collectors  and  by  contact 
resistance  between  the  cell  components.  Ohmic  losses  can  be  cal¬ 
culated  as  follows: 


Elohm  —  Rtoti  (66) 


where  Rtot  is  the  total  area  specific  resistance  and  can  be  calculated 
as  follows  [142]: 


where  <5  is  the  thickness,  and  a,  c  and  e  denote  the  anode,  cathode 
and  electrolyte,  respectively,  and  a  is  the  conductivity  of  electrodes 
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and  electrolyte  which  is  defined  as  a  function  of  temperature  as 
follows  [143]: 


ae  =  3.34  x  104 


10300\ 

) 


a  a  = 


95  x  106 


exp  - 


1150  \ 
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(7C  = 


f  42  x  10 

' 


/  1200\ 

6XP(-  —  J 


(68) 


The  above  equations  are  valid  for  the  YSZ  electrolyte  material 
[144]. 

Hajimolana  and  Soroush  [22]  and  Qi  et  al.  [34,43]  assumed 
a  constant  value  for  specific  resistance  while  Costamagna  et  al. 
[26],  Stiller  et  al.  [28]  and  other  many  research  studies  calcu¬ 
lated  the  anode/electrolyte/cathode  resistance  [25,26,28,69,110, 
133,145]. 

In  order  to  evaluate  ohmic  losses,  some  authors  [28,59]  used  the 
Nisancioglu’s  expression,  Eqs.  (69)-(74)  [146],  with  a  correction  for 
the  ionic/electric  resistivity  dependence  on  temperature,  Eq.  (74) 
giving  by: 


A  i/  jjAsliceReq,j 

vohm,i  —  t  ? 

^slice 

(69) 

where 

Req  =  Req  \  +  Req, 2? 

(70) 

Req,  1  — 


((Pa/^a)  +(Pc/^c)  cosh(Ze)  +  papc/ <Wc(2  +  Ze  sinh(Ze))) 


2(1  /pe&e)V2{Pa/Za  +  Pc/&cf/2  sinh(Ze) 


(71) 


Req,2  — 


\J  Picm &icm  ( Pc  / &c ) 


2  tanh(LjCm/2y/ (Pc/^c)/ Picm&icm) 


(72) 

(73) 


Pi  = 


(74) 


In  these  equations  A  is  the  active  area,  jj  the  current  density, 
Reqj  is  equivalent  ohmic  resistance,  L  is  the  slice  length,  p  ohmic 
resistance,  <5  thickness  and  C  resistivity. 

Several  SOFC  models  have  been  developed  to  study  polarization 
losses  within  the  PEN  structure. 

Hirano  et  al.  [146]  presented  a  model  in  order  to  investigate 
the  effect  of  concentration/ohmic  polarization  on  electrochemical 
model  by  assuming  that  the  activation  polarization  is  negligible. 

Burt  et  al.  [27],  Ji  et  al.  [29],  Chan  et  al.  [30],  Izzo  Jr 
et  al.  [38]  and  other  many  papers  like  [25,27,29,30,38,41,66, 
74,77,78,83,97,108,110,115,117,120,134,147-151]  performed 
the  three  activatin/concentration/ohmic  polarizations  in  their 
model. 

Tanaka  et  al.  [152],  Aguiar  et  al.  [108]  and  Achenbach  [116] 
considered  only  the  activation  polarization  in  their  model. 

In  the  Aguiar  et  al.  [108]  survey,  it  was  assumed  that  except 
for  the  shift  reaction  the  principle  electrochemical  reactions  are 
considered  kinetically  controlled,  and  therefore  only  the  activation 
overpotentials  were  taken  into  account.  It  was  assumed  that  the 
other  overpotentials  are  small  and  can  be  ignored. 

Chaisantikulwat  et  al.  [39],  Ho  et  al.  [45],  Bhattacharyya  et  al. 
[54],  Janardhanan  et  al.  [67],  Danilov  et  al.  [99]  and  Wang  et  al. 


[98]  accounted  both  activation  and  ohmic  losses  in  their  model  and 
neglected  the  concentration  polarization. 

Activatioin  and  ohmic  polarization  were  accounted  by  Ota  et  al. 
[34].  In  the  study  it  was  also  assumed  that  concentarion  gradient 
in  the  gas  phase  exists  only  in  the  porous  electrodes. 

PadulleTs  et  al.  [153]  neglected  the  activation/concentration 
polarization  and  only  considered  to  ohmic  loss  in  their  model  while 
Yakabe  et  al.  [42]  considered  only  concentration  polarization  in 
their  model. 

In  recent  years,  several  SOFC  models  have  been  developed  by 
Chan  et  al.  [30],  Sanchez  et  al.  [59],  Nam  et  al.  [64],  Klein  et  al.  [90] 
and  other  studies  in  [30,55,59,62-64,85,86,90,136,148,150,154]  to 
investigate  reactions  and  transport  phenomena  within  the  PEN 
structure.  These  models  differ  widely  in  terms  of  their  complexity 
and  comprehensiveness. 

Chan  et  al.  [30]  performed  a  complete  polarization  model.  In 
this  paper  the  effect  of  electrodes  and  electrolyte  thicknesses  was 
studied.  It  was  observed  that  the  sensitivity  of  cell  voltage  due  to 
change  of  electrolyte  thickness  is  the  highest,  then  sensitivity  of  cell 
voltage  due  to  the  change  of  cathode  thickness  and  finally  the  sensi¬ 
tivity  of  cell  voltage  due  to  the  change  of  anode  thickness.  Reducing 
all  three  component  thicknesses  causes,  in  general,  reduced  sensi¬ 
tivity  strength  with  an  improved  operating  current  density  range. 
However;  considering  both  current  density  and  PEN  thicknesses 
is  significant  and  needs  an  optimization  between  both  factors  in 
order  to  obtain  a  high  performance  while  in  this  paper  both  fac¬ 
tors  was  issued.  A  cathode-supported  cell  is  no  way  better  than  an 
anode-supported  cell  for  high  performance  fuel  cell  design  even 
under  elevated  pressure  conditions  at  the  cathode  compartment 
to  compensate  for  the  serious  cathode  concentration  overpotential 
loss.  It  was  illustrated  that  the  cell  potential  losses  is  mostly  due  to 
anode  side  of  the  cell  in  an  anode-supported  SOFC. 

Nam  et  al.  [64]  illustrated  a  comprehensive  micro-scale  model 
for  explaining  the  transport  phenomena  and  reaction  in  SOFCs. 
They  showed  that  the  potential  loss  in  cathode  was  considerable  at 
all  current  densities  while  the  potential  loss  in  anode  was  negligible 
at  low  current  densities  and  important  at  higher  current  densities. 
Based  on  their  research,  it  was  observed  that  the  size  of  particle 
has  the  most  significant  impact  on  the  PEN  performance  optimiza¬ 
tion.  By  increasing  the  length  of  TPB  in  the  electrodes  and  smaller 
particle  size,  the  activation  overpotential  will  be  reduced,  while  at 
the  same  time  by  decreasing  Knudsen  diffusivity  and  flow  perme¬ 
ability  the  mass  transport  resistance  increases.  On  the  other  hand, 
the  volume  fraction  of  electronic  phase  of  around  0.45  was  found 
optimal  for  both  anode  and  cathode  for  the  same  ionic  and  elec¬ 
tronic  particle  size.  Higher  porosity  in  anode  and  lower  porosity 
in  cathode  gives  a  better  PEN  performance.  The  porosity,  particle 
size,  and  thickness  of  electrodes  were  found  to  have  conflicting 
effects  on  the  PEN  performance.  For  example,  by  decreasing  the 
particle  size  of  electrodes  the  efficiency  of  electrochemical  reac¬ 
tion  increases  while  the  efficiency  of  mass  transport  decreases. 
An  optimal  set  of  micro-structural  parameters,  which  resulted  in 
best  PEN  performance,  will  not  automatically  result  in  best  SOFC 
performance.  Micro-structural  optimization  of  SOFC  performance 
requires  consideration  on  the  interconnect  rib  geometry  and  the 
flow  field  effects,  in  addition  to  the  detailed  micro-scale  calcula¬ 
tion  in  PEN.  In  the  paper  also  the  effect  of  electrodes  thickness  was 
investigated  in  details.  It  was  also  depicted  that  with  a  larger  of 
electrode  thickness  the  performance  increases,  however;  at  higher 
current  density  the  effect  of  electrode  thickness  will  be  negligi¬ 
ble. 

Jiang  et  al.  [148]  modeled  a  one-dimensional  dynamic  model 
of  a  TSOFC  stack.  It  was  found  that  higher  operating  temperature 
decreases  both  the  Nernst  potential  and  the  irreversible  losses, 
resulting  in  an  initial  increase  then  a  decrease  in  cell  efficiency. 
In  this  paper  the  resistivity  and  ohmic  losses  in  different  layers  of 
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tubular  SOFC  were  calculated  as  follows: 

p  =  0.00298  e~1392/r  fuel  electrode 
p  =  0.00294  e10350/7  electrolyte 
p  =  0.0081 14  e600/7  air  electrode 
p  =  0.12568  e4690/7  interconnect 

R  =  pLA 

Jiang  et  al.  also  studied  the  effect  of  fuel  cell  size  like  fuel  cell 
diameter  and  length.  It  was  represented  that  by  an  increase  of  cell 
diameter  the  power  increase  due  to  a  larger  activation  area  at  the 
same  time  and  also  the  ohmic  loss  increase  due  to  longer  current 
path  length. 

Ni  et  al.  [66]  developed  an  electrochemical  model  for  studying 
the  ammonia  (NH3)-fed  solid  oxide  fuel  cells  with  proton¬ 
conducting  electrolyte  (SOFC-H)  and  oxygen  ion-conducting 
electrolyte  (SOFC-O).  It  was  observed  that  the  actual  performance 
of  the  NH3-fed  SOFC-H  is  significantly  lower  than  SOFC-O,  mainly 
due  to  higher  ohmic  overpotential  of  the  SOFC-H  electrolyte.  It  was 
also  presented  that  compared  with  the  NH3-fed  SOFC-H,  the  SOFC- 
O  has  higher  anode  concentration  overpotential  and  lower  cathode 
concentration  overpotential.  The  effects  of  temperature  and  elec¬ 
trode  porosity  on  concentration  overpotentials  were  also  studied 
in  order  to  identify  possible  methods  for  improvement  of  SOFC 
performance.  This  study  observed  that  the  use  of  different  elec¬ 
trolytes  not  only  causes  different  ion  conduction  characteristics  at 
the  electrolyte,  but  also  significantly  influences  the  concentration 
overpotentials  at  the  electrodes. 

Ni  et  al.  [41  ]  investigated  the  effect  of  micro-structures  of  elec¬ 
trodes  on  SOFC  performance.  A  mathematical  model  for  a  planar 
SOFC  was  developed  for  studying  the  effects  of  varies  parameters  on 
the  fuel  cell  performance.  Nam  et  al.  [64]  work,  they  revealed  that 
both  porosity  grading  and  pore  size  grading  can  effectively  improve 
the  SOFC  performance.  Based  on  this  survey,  at  low  current  density 
SOFC  potential  decreases  with  increasing  porosity,  while  at  high 
current  density,  optimal  porosity  was  found  and  the  optimal  poros¬ 
ity  increases  with  increasing  current  density.  Moreover;  the  effect 
of  pore  size  and  porosity  of  electrodes  on  the  overpotentials  was 
studied  while  this  effect  was  not  considered  by  Nam  et  al.  [64]  work. 
Based  on  their  results,  variation  of  porosity  has  an  inverse  influence 
on  activation  and  concentration  overpotentials.  Similar  to  poros¬ 
ity,  the  electrode  pore  size  has  opposite  effects  on  activation  and 
concentration  overpotentials.  The  mechanisms  of  coupled  trans¬ 
port  and  chemical  reactions  in  the  porous  electrodes  and  evaluated 
possible  ways  to  improve  SOFC  performance  were  also  studied.  The 
model  developed  in  the  paper  can  be  effectively  applied  to  the  SOFC 
design  optimization  to  maximize  the  energy  efficiency  of  CH4  fed 
SOFC  for  clean  power  generation.  In  addition;  the  impact  of  CH4 
fed  SOFC  on  overpotentials  was  studied.  Regarding  this  study,  for 
SOFC  fed  by  CH4,  all  the  overpotentials  decrease  with  increasing 
temperature  while  in  H2  fed  SOFC  the  concentration  overpoten¬ 
tial  slightly  increases  with  increasing  temperature.  The  effect  of 
methane  reforming  and  water-gas  shift  reactions  is  significant  at 
high  temperature,  resulting  in  high  rate  of  H2  production  and  high 
molar  ratio  of  H2  and  H20. 

Arpino  et  al.  [74]  investigated  the  effect  of  current  density  and 
fuel  utilization  on  the  overpotential  losses.  Based  on  their  research, 
both  concentration  and  activation  polarizations  at  the  anodic  com¬ 
partment  increase  significantly  for  a  high  value  of  current  density. 
It  was  also  observed  that  there  are  no  important  gradients  in  the 
anodic  concentration  polarization  at  lower  current  density  as  the 
inlet  hydrogen  velocity  varies. 

Table  6  illustrates  some  papers  that  concerned  on  PEN  in  their 
model. 


(75) 

(76) 


3.4.3.  Mass  transfer  model 

At  the  electrodes,  hydrogen  and  oxygen  in  the  case  of  pure 
hydrogen  as  a  fuel,  diffuse  from  anode  and  cathode  sides,  respec¬ 
tively,  into  the  triple  phase  boundary  for  occurring  electrochemical 
reaction,  and  then  after  the  products  diffuse  into  the  anode  side.  The 
diffusion  of  components  from  fuel  channel  into  the  TPB  depends  on 
the  type  of  the  fuel  fed  to  the  SOFCs.  Mass  transport  models  inside 
the  electrodes  are  giving  an  accurate  prediction  for  gas  concen¬ 
tration  at  the  anode-electrolyte-cathode  interface.  Therefore;  an 
accurate  modeling  of  this  phenomenon  is  imperative  for  develop¬ 
ment  of  better  fuel  cell  designs. 

A  mass  balance  on  oxygen  inside  the  cathode-side  diffusion 
layer  yields: 


A  cat  9 


vn0  — 


R 


o2 


Ac 


(77) 


Acat  is  the  thickness  of  the  cathode  layer,  R  is  the  universal  gas 
constant,  P JPB  is  the  pressure  of  component  at  TPB,  Ac  is  the  surface 
area  of  cathode  side  and  Tctrpe  is  the  cell  tube  temperature  at  TPB. 

By  assuming  that  synthesis  gas  is  as  a  fuel  fed  to  SOFC,  mass  bal¬ 
ances  on  hydrogen  and  water  vapor  inside  the  anode-side  diffusion 
layer  yield  are  given  as  follows: 


Aano  9 
—R~di 
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Aano  is  the  thickness  of  the  anode  layer,  Aa  is  the  surface  area  of 
anode  side. 


3.4.4.  Heat  transfer  model 

Cell  temperatures  influence  the  electrochemical  model  and 
these  temperatures  affect  the  local  driving  voltage,  polarizations, 
and  heat  generation  within  the  cell. 

(i)  An  energy  balance  inside  the  SS4  for  TSOFC  leads  to: 

mct  at  Ct  =  ^°nv  +  ®c°nd  +  Qrad  +  Q diff  +  ®elc  (80) 

Convection  heat  transfer  inside  the  SS4  for  TSOFC  is  between 
inside  the  wall  of  the  cell  tube  and  air  flow  inside  the  SS2  and 
outside  the  wall  of  cell  tube  and  fuel  flow  which  is  given  by: 

QcCL  =  ha, (JZ  ~  Tcti )  +  hct0(Tp  -  Tct0)  (81) 

Thermal  radiation  within  tubular  SOFC  in  SS4  occurs  between 
injection  tube  and  inside  the  cell  tube  walls  which  is  given  by: 

Va,  =  if-K  -  TcV  (82) 

lxrad 

For  both  TSOFC  and  PSOFC  in  SS4,  02  and  H2  diffuse  into 
the  cathode  and  anode,  respectively,  in  order  to  participate  in 
electrochemical  reaction,  and  H20  is  the  electrochemical  pro¬ 
duction  that  diffuses  from  cell  tube  into  the  fuel  channel.  Thus 
the  diffusional  heat  transfer  leads  to: 

os§='f2VNiHi  (83) 

In  the  case  of  synthesis  gas  as  a  fuel,  02,  H2,  and  CO  diffuse 
into  the  cell  tube,  and  H20  diffuses  into  the  fuel  channel. 

In  the  case  of  electrochemical  heating  (Qe/ec),  heat  gener¬ 
ates  by  the  reversible  and  irreversible  processes  inside  the 
electrodes  and  electrolyte  within  the  triple  phase  boundary. 
Reversible  heat  generation  is  due  to  entropy  transfer  by  the 
electrochemical  reactions  at  the  cathode  and  anode  side  which 
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Table  6 

Selected  studies  concerning  on  PEN  model  of  SOFCs. 


Author 


Cell  type 


Objective 


Method 


Short  coming 


Assumption 


F.  Zhao  etal.  [158]  PSOFC 


N.  Autissier  et  al.  [35]  PSOFC 


Shi  etal.  [86]  PSOFC 


Xie  and  Xue.  [87]  anode-supported  solid 

oxide  fuel  cell 


To  anticipate  the 
concentration  polarization 
in  detail. 


To  predict  current  density, 
flow,  temperature  and 
concentration  fields  in 
order  to  compare  and 
optimize  repeat  element 
geometry  for  a  whole  stack. 


To  investigate  the  intricate 
interdependency  among 
the  ionic  conduction, 
electronic  conduction, 
multi-component  species 
transport,  electrochemical 
reaction  processes  and 
electrode  microstructure 
for  intermediate 
temperatures  operation, 
between  750  and  850  °C 


To  de-convolute  the 
complicated  transient 
interactions  among 
operating  conditions, 
internal  multi-physics 
processes,  and  SOFC 
current  responses. 


Porous  media  gas  phase 
transport  processes  based 
on  Fields  model  was 
considered  for  predicting 
the  concentration  over 
potentials. 

A  commercial  CFD  tool  was 
used,  solving  mass, 
momentum  and  energy 
equations;  whereas 
chemical  kinetic  equations 
were  computed  from 
external  sub-routines, 
Radiative  heat  transfer  was 
accounted  between  inner 
surfaces. 

The  model  was  calibrated 
using  experimental 
polarization  curves  and 
then  validated  by 
comparing  each  cell 
component  polarizations 
(anodic,  cathodic  and 
electrolyte) 


The  model  was  tacked  into 
account  the  complex 
transport  phenomena 
within  PEN  assembly, 
including  charge 
(ion/electron)  migration 
and  species  transport. 


The  models  are  limited  to 
particular  cells  and  cannot 
describe  the  effects  of 
operating  and  design 
parameters  in  detail. 

The  model  parameters  are 
not  determined  from 
experiments  or  from 
literature,  they  were  only 
estimated. 


-Convection  flux  was 
neglected  in  the  porous 
electrode  compared  to 
diffusion. 


-  Pressure  gradients  in  the 
porous  electrode  were 
neglected. 

-  All  physical  properties 
were  evaluated  at  a  fixed 
cell  temperature. 

The  pressure  gradient 
within  porous  electrodes 
was  neglected 


ID 


ID  steady  state, 
laminar  for  both  fuel 
and  air  flow 


2D  steady  state, 
isothermal 


Isothermal  2-D 
transient 


is  exothermic  and  slightly  endothermic,  respectively  [155] 
which  is  given  by: 


Qrev,a  —  TctTPB{2Sy[2Q 
Qrev,c  —  TctTpB(2So2- 


2Sh2  +4Se-)4jF 

(84) 

2Sq2  4Se-)4JF 

(85) 

where  Se  is  the  molar  entropies  of  the  products  and  reactants 
at  the  operating  temperature  and  pressure. 

Irreversible  heat  generation  at  the  anode  and  electrolyte 
region  is  caused  by  anodic  and  cathodic  overpotentials  which 
are  given  by: 


All  the  heat  sources  generated  by  electrochemical  reac¬ 
tion  are  determined  by  Bove  et  al.  [32],  Klein  et  al. 
[90],  Ni  et  al.  [103],  Young  et  al.  [104]  and  studies  in 
[25,32,39,90,103,104,106,126,129,131,157,158]. 

(ii)  An  energy  balance  inside  the  SS4  for  PSOFC  leads  to: 

mct  Ct  =  ^°nv  +  Q°nd  +  +  +  ^lc  (89) 

Convection  heat  transfer  inside  the  SS4  for  PSOFC  is  between 
the  wall  of  the  cathode  side  and  air  flow  inside  the  SSI  and  the 
wall  of  the  anode  side  and  fuel  flow  which  is  given  by: 

Qconv  =  het, c0 ( Tair  ~  Tct,c0)  +  ^ct, a0 ( ^fUel  ~  ^cTao)  (90) 


Qirr,a  —j^la 
Qirr,c  =Jr/c 


(86) 

(87) 


Heat  also  generates  by  the  ohmic  heating  in  the  electrodes 
and  electrolyte  due  to  existing  of  resistances  to  flow  of  ions 
and  electrons  which  is  given  by: 


Thermal  radiation  within  planar  SOFC  inside  the  SS4  occurs 
between  cathode  side  of  cell  tube  and  air  side  interconnector 
walls,  and  also  between  anode  side  of  the  cell  tube  and  fuel  side 
interconnector  walls  which  is  given  by: 


_ (rrano4  _ rr4  n  .  _ rrrcat4  _ y4  n 

d  ,^Jint ,i  ict,a0^i_  d  Ict,c0l 

lxrad  iyrad 


(91) 


Qohm  =  Pei2  (88) 

where  pe  is  the  ionic  resistivity  of  the  electrolyte.  It  is  reported 
that  ohmic  heating  is  about  2.37-4.1%  of  the  total  heat  pro¬ 
duced  in  an  electrolyte-supported  planar  SOFC  [156].  However, 
Andreassi  et  al.  [129]  assumed  that  the  ohmic  heating  in  the 
porous  electrodes  is  negligible  because  of  high  electrical  con¬ 
ductivity  as  compared  to  the  ionic  conductivity. 


In  some  works,  radiation  heat  transfer  was  neglected  in  the 
planar  SOFC  because  in  these  works  it  was  shown  that  this  type 
of  heat  transfer  has  only  a  negligible  effect  on  the  temperature 
field  within  the  electrode  and  electrolyte  [155,159,160]. 

Large  variations  in  the  cell  performance  occurs  when  only 
conductive  and  convection  heat  transfer  among  the  fuel  cell 
stack  is  considered,  while  by  adding  the  radiation  heat  to  the 
model,  the  uniformity  of  temperature  promotes  within  the 
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Table  7 

Selected  studies  concerning  on  novel  solid  oxide  fuel  models. 


1909 


Author 

Objective 

Method 

Results  and  advantages  of  the  model 

Entchev,  E  et  al.  [192] 

To  predict  SOFC  performance  while 
supplying  both  heat  and  power  to  a 
residence. 

Applied  adaptive  neuro-fuzzy 
inference  system  (ANFIS) 
techniques  and  artificial  neural 
network  (ANN) 

-The  prediction  of  both  ANN  and  ANFIS  models 
was  well  in  agreement  with  variety  of 
experimental  data. 

-Compared  to  the  conventional  models  the 

ANN  and  ANFIS  models  are  able  to  predict  and 
optimize  system  performance  faster  and 
deliver  better  results  in  many  instances. 

-the  ANFIS  model  predicts  the  fuel  cell  current 
better  than  ANN  model. 

-The  ANN  model  predict  the  fuel  cell  voltage 
better  than  the  ANFIS  model. 

Wu  et  al.  [197] 

To  describe  the  nonlinear 
temperature  and  voltage  dynamic 
properties  of  the  SOFC  system 

It  was  applied  Takagi-Sugeno 
(T-S)  fuzzy  model. 

-  The  relationship  between  inputs  and  outputs 
of  the  SOFC  system  can  be  predicted  while  the 
complicated  inner  structure  was  ignored. 

Arriagada  et  al.  [198] 

To  evaluating  the  SOFC 
performance 

Used  an  Artificial  neural  network 
method 

-The  SOFC  model  based  on  ANN  was  well  in 
agreement  with  the  physical  model. 

-ANN  model  was  much  easier  and  faster  to  use 
compared  to  mathematical  model. 

Yang  et  al.  [199] 

To  control  fuel  cell  temperature 

A  modified  Takagi-Sugeno  (T-S) 
fuzzy  modeling  study  was 
presented. 

-The  modified  T-S  fuzzy  model  simplified  the 
complicated  mathematical  model  and  quickly 
made  an  accuracy  relationship  between  inputs 
and  outputs  of  the  system. 

-Based  on  the  results  the  modified  T-S  fuzzy 
yielded  a  high  accuracy  over  a  relatively  wide 
operating  temperature  range 

Wu  et  al.  [200] 

To  facilitate  a  valid  control  strategy 
design 

A  radial  basis  function  (RBF) 
neural  network  based  on  a 
genetic  algorithm  (GA)  was 
applied. 

-The  results  were  congruence  in  agreement 
with  simulation  results. 

-  It  was  possible  to  design  an  online  controller 
of  a  SOFC  stack  based  on  this  GA-RBF  neural 
network  identification  model. 

Hua  et  al.  [201  ] 

-To  facilitate  a  valid  control 
strategy  design. 

-  Trying  to  avoid  the  complexities 
of  model  based  on  conversion  laws 

Performed  a  black-box  model  of 
the  SOFC  based  on  least  squares 
support  vector  machine 
(LS-SVM) 

-The  LS-SVM  was  superior  to  the  conventional 
RBFNN  in  predicting  stack  voltage  with 
different  fuel  utilizations. 

stack,  thus  leads  to  more  uniform  the  cell  voltage.  The  varia¬ 
tion  of  cell-to-cell  voltage  is  less  than  0.2%  in  the  model  that  is 
involved  radiation  heat  transfer  and  0.3%  in  the  model  without 
radiation  [27].  In  the  case  of  indirect  internal  reforming,  that 
radiative  heat  transfer  accounts  for  up  to  79%  of  the  total  heat 
transfer  between  the  solid  structure  and  the  reformer  [108]. 

3.4.5.  Reforming  and  shift  reactions 

Reforming  reaction  of  hydrocarbons  such  as  methane  or  natu¬ 
ral  gas  is  a  process  for  producing  the  hydrogen  required  by  fuel 
cell.  This  reaction  is  predominately  performed  between  600  and 
900  °C  and  is  thus  compatible  with  SOFCs.  The  heat  produced  in  a 
fuel  cell  can  provide  the  heat  for  the  endothermic  reforming  reac¬ 
tion.  Therefore,  with  elaborating  an  appropriate  catalyst  using  in 
anode  manufacturing,  internal  reforming  reaction  in  fuel  cell  can 
eliminate  the  requirement  for  indirect  internal  fuel  reformer  and  is 
expected  to  simplify  the  overall  system  design  and  finally  the  eco¬ 
nomical  cost  of  these  power  systems  will  decrease  substantially.  By 
manipulating  the  internal  reforming  reaction  the  fuel  cell  tempera¬ 
ture  can  decrease  to  a  suitable  temperature.  There  are  two  alternate 
internal  reforming  within  a  SOFC:  indirect  (IIR-SOFC)  and  direct 
(DIR-SOFC)  internal  reforming.  In  the  former  process  the  reformer 
part  is  separate  but  adjacent  to  the  anode  cell.  In  the  second  process, 
methane  is  fed  directly  into  the  cell  and  the  reaction  occurs  on  the 
anode.  However,  this  process  has  two  major  problems.  One  prob¬ 
lem  is  that  internal  reforming  is  highly  endothermic  and  thus  cell 


temperature  drops  initially  and  then  increases  again,  hence  gener¬ 
ate  large  temperature  gradients  across  the  cell,  as  well  as  a  drop  in 
cell  performance.  Deactivation  of  the  anode  by  carbon  formation  at 
intermediate  temperature  and  high  pressure  is  the  second  problem 
of  the  DIR-SOFC  which  reduces  the  system  efficiency.  The  problem, 
however,  can  be  solved  by  adding  H20  or  C02  into  the  syngas  fuel. 
The  complete  steam  reforming  reaction  takes  place  according  to 
Eq.  (92) 

CH4  +  2H20  -*  4H2  +  C02  (92) 

Reforming  reaction  mainly  occurs  in  the  diffusional  layer  of 
anode  or  in  the  electrochemically  active  layer.  But  the  rate  of  the 
reaction  in  the  latter  layer  is  assumed  to  be  much  lower  than  the  for¬ 
mer  layer  because  of  competition  with  the  electrochemical  reaction 
which  includes  much  more  hydrogen. 

Eq.  (92)  is  the  combination  of  two  independent  reactions: 
methane  steam  reforming,  Eq.  (93),  and  water  gas  shifting,  Eq.  (96), 
which  are  highly  endothermic  and  slightly  exothermic,  respec¬ 
tively: 

CH4  +  H20^3H2  +  C0  (93) 

CO  +  H20^  H2  +  C02  (94) 

The  shift  reaction  assumed  to  tacke  place  always  at  the  anode 
layer. 
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Although  Eqs.  (93)  and  (94)  occur,  the  following  carbon  forming 
reactions  may  also  appear  in  the  anode  and  the  pre  reformer: 

CH4^C  +  2H2  (95) 

2CO  ->  C  +  C02  (96) 

These  two  reactions  are  endothermic  and  exothermic,  respec¬ 
tively.  If  any  of  these  two  reactions  take  place  at  the  anode, 
deposition  of  carbon  particles  over  the  anode  surface  will  deac¬ 
tivate  the  catalyst  and  thus  reduce  the  cell  performance.  At  high 
temperature  around  1000K,  methane  is  more  likely  to  react  with 
steam  to  form  hydrogen  according  to  Eq.  (92)  than  to  dissociate 
into  atomic  carbon  and  water  [150].  This  reaction  increases  when 
the  steam  concentration  is  high.  For  carbon  monoxide,  decomposi¬ 
tion  occurs  when  the  steam  concentration  is  not  high  enough.  Shift 
reaction  is  clearly  dominant  by  the  excess  steam  in  the  fuel  stream. 

Many  researchers  like  Hajimolana  and  Soroush  [22], 
Ni  et  al.  [41],  Qi  et  al.  [43],  Nikooyeh  et  al.  [48],  Suwan- 
warangkul  et  al.  [52]  and  also  other  researchers  in 
[22,41,43,48,52,65,67,69,91,97,98,110,150,161,162]  considered 
DIR-SOFC  in  their  model. 

Yingwei  et  al.  [97]  observed  that  direct  internal  reforming  pro¬ 
cess  has  significant  effect  on  the  distributions  of  the  anode  gas 
composition,  cell  temperature,  solid  temperature  gradient,  and  cur¬ 
rent  density  along  the  cell  length  direction. 

Dokamaingam  et  al.  [77,78],  Brus  et  al.  [121],  Aguiar  et  al. 
[108,109]  modeled  IIR-SOFC  in  their  system. 


3.4.5 A.  Equilibrium  model.  Few  authors  considered  that  the 
reforming  and  shift  reaction  are  equilibrium  inside  the  cell 
[163,164].  The  constant  equilibrium  for  reforming  and  shift 
reactions  yield  with  the  [99,100]  equations  which  depend  on  tem¬ 
perature  exponentially  as  depicted  in  Eq.  (100),  where  subscript  i 
stands  for  both  reforming  and  shift  reactions.  The  constant  values 
are  reported  in  [163,164]  which  are  given  by: 


Keqjef  : 
^eq,  shift 


Ph/co 

/Wh2o 

Ph2Pco2 

Ph2oPco 


(97) 


log  KeqJ  —  aiPct  +  bi^ct  +  Ci^ct  +  d{Tct  +  6/ 


(98) 


which  both  the  reforming  and  shift  reactions  occur  at  the  surface 
of  the  anode  side. 

Suwanwarangkul  et  al.  [61]  calculated  the  shift  reaction  equi¬ 
librium  constant  as  follows: 


I<eq, shift  =  exp(-0.2935 Z3  +  0.6351Z2  +  4.1788Z  +  0.3169)  (99) 


where 


Z  = 


1000 

W) 


(100) 


However,  equilibrium  model  is  not  a  common  assumption  for 
refroming  reaction  in  the  papers. 


3.4.52.  Kinetic  model.  Some  researchers  claim  that  equilibrium 
cannot  be  assumed  for  the  reforming  reaction  inside  the  solid  oxide 
fuel  cell  due  to  a  slow  reforming  kinetics  compared  with  both 
shifting  and  hydrogen  electrochemical  oxidation.  Therefore,  the 
following  general  expression  is  used  to  evaluate  the  reforming  rate 
and  depends  on  reactants  concentration  and  temperature. 

rCH4  =  •  Ph20  eXP  )  (10^) 

where  Ea  is  the  activation  energy  and  X  is  the  pre-exponential  fac¬ 
tor. 


Hajimolana  and  Soroush  [22],  Qi  et  al.  [43],  Nehter  et  al.  [75] 
and  Aguiar  et  al.  [108]  assumed  a  kinetic  reforming  reaction  in  their 
model. 

3.5.  Model  equations  for  subsystem  5  (SS5) 

This  subsystem  for  TSOFC  is  the  fuel  that  occupies  and  flows  in 
the  space  on  the  anode  side  of  the  cell  tube  and  for  PSOFC  is  the 
fuel  channel. 


3.5 A.  Mass  transfer  model 

The  multi-component  mass  transport  in  SOFC  anodes  is  mainly 
coupled  with  the  bulk  chemical  reactions,  diffusion  and  convection 
of  species  while  at  the  cathode  side  there  is  no  chemical  reaction. 
A  total  mass  balance  on  the  fuel  and  component  mass  balances  on 
the  fuel  components  for  both  TSOFC  and  PSOFC  leads  to: 


(102) 

(103) 


Wj  is  the  rate  of  producing  or  consuming  of  components. 


3.5.2.  Heat  transfer  model 

(i)  An  energy  balance  for  the  fuel  gas  inside  SS5  for  TSOFC  results 
in: 


dp 


anoijano 

fuelnfuel 


dt 


VAp 


anoi.anoj^ano  ,  rtano 


'fuelnfuel1 


fuel 


+  dC 


(104) 


In  SS5,  the  convection  heat  transfer  for  TSOFC  and  PSOFC 
occurs  between  outside  of  the  cell  tube  wall  and  the  fuel  flow 
inside  the  SS5  which  is  given  by: 


QcaZ  =  hct0(Tct0-Tp)  (105) 

Qcmo  jn  tj1js  subsystem  for  both  TSOFC  and  PSOFC  includes 
radiation  heat  transfer  between  walls  and  gas  components. 
However,  in  most  of  the  works,  this  type  of  radiation  was 
neglected  due  to  the  lack  of  sufficiently  accurate  properties  of 
gases. 

Fuel  components  including  H20,  CH4,  CO,  C02  and  possibly 
other  hydrocarbons  are  very  strange  spectral,  temperature  and 
pressure  dependence  and  so  can  be  treated  as  transparent  [93]. 

This  heat  flow  per  wall  unit  area  is  evaluated  by  using  Eq. 
(106)  where  the  emissivity  of  the  mixture  of  gases  sg  is  cal¬ 
culated  as  proposed  by  Hadvig  in  [165],  Eq.  (107)  and  its 
absorptivity  ag  as  shown  in  Eq.  (110)  [166].  Subscripts  w  and 
g  stands  are  for  wall  and  gas  in  Eqs.  (106)-(108)  and  i  for  the 
slice  being  analyzed. 


<Jrad,i  Ugyy  _  sw)a^£gTS>  agTwf 

Sg  =  sC02  -  £H20  +  £C02£H20 

/  t  \  0.65  /  T  \ 

ag  =  eC02  (  A)  +  eH2o(^)  -eC02eH20 


(106) 

(107) 

(108) 


Few  studies  considered  radiation  heat  transfer  between  walls 
and  gas  components  in  their  model  [59]. 

(ii)  An  energy  balance  for  the  fuel  gas  inside  SS5  for  PSOFC  results 
in: 


dp 


anoijano 

fuelnfuel 


dt 


V7  a  r.ano1  ,ano  jrano  ,  r\ano 
V  APfuelUfuelHfuel  +  ^ 


(109) 
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Fig.  5.  The  behavior  of  Fuel  Cell  Voltage  at  steady-state  and  transient  condition. 


In  SS5,  the  convection  heat  transfer  for  PSOFC  occurs  between 
both  the  anode  side  of  the  cell  tube  wall  and  fuel  side  intercon¬ 
nector  and  the  fuel  flow  inside  the  SS5  which  is  given  by: 


that  SOFC  can  use  many  common  hydrocarbon  fuels  such  as  natural 
gas,  diesel,  gasoline  and  alcohol  without  the  need  to  reform  the  fuel 
into  pure  hydrogen.  In  other  fuel  cells,  such  as  the  polymer  elec¬ 
trolyte  fuel  cell,  which  are  fueled  with  pure  hydrogen,  the  carbon 
dioxide  is  a  poison. 

Hydrogen  was  considered  as  a  fuel  for  SOFCs  by  Julian  et  al.  [7], 
Wang  et  al.  [126],  Yuan  et  al.  [168]. 

Hajimolana  and  Soroush  [22],  Autissier  et  al.  [35],  Aloui  et  al. 
[69],  Aguiar  et  al.  [108],  Cui  et  al.  [101],  Janardhanan  et  al.  [157], 
Fryda  et  al.  [169]  and  Colpan  et  al.  [170]  considered  hydrocarbon 
fuel  fed  to  SOFCs. 

Many  authors  like  Suwanwarangkul  et  al.  [61  ],  Aloui  et  al.  [69], 
Xie  et  al.  [86]  and  other  researchers  in  [61,69,86,171-179]  tacked 
into  account  biomass  as  a  fuel  for  SOFCs  because  it  is  abundant  and 
renewable. 

Hydrogen  sulfide  (H2S)  is  a  byproduct  of  natural  gas  which  can 
also  use  as  a  fuel  in  the  solid  oxide  fuel  cells.  There  are  two  possible 
overall  reactions  that  use  H2S  directly  in  a  SOFC: 

H2S+1o2^H20+1s2  (112) 

H2S+^02^H20  +  S02  (113) 


Qconv  —  hct,a0(Tct,a0 


jch,a\  ,  l  fj, 
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Tch,as 
fuel  j 


(110) 


The  temperature  distribution  within  the  fuel  cell  is  vital  to 
SOFC  operation,  affecting  the  reaction  kinetics  and  electro¬ 
chemical  model  and  these  temperatures  influence  the  voltage, 
polarization,  heat  generation  within  the  cell,  thus  the  over¬ 
all  cell  performance.  Based  on  the  literatures,  when  only 
convection/conduction  or  convection/radiation  heat  transfer 
is  considered  large  variations  in  cell  temperature  distribu¬ 
tion  and  then  in  the  cell  performance  was  observed.  In 
Fig.  6,  there  is  a  comparison  between  three  different  mod¬ 
els  regarding  to  distribution  of  cell  temperature  within  the 
cell  length.  A  complete  heating  components  including  radi¬ 
ation/convection/conduction/electrochemical  and  diffusional 
heating  is  involved  in  the  T.  Ota  et  al.  [167]  model,  while  R. 
Suwanwarangkul  et  al.  [61  ]  studied  two  different  energy  mod¬ 
els,  first  a  model  including  only  conduction,  convection  and 
radiation  and  in  the  second  model  the  radiation  heat  transfer  is 
ignored.  It  is  observed  that  the  temperature  difference  between 
the  minimum  temperature  and  maximum  temperature  within 
the  cell  length  is  87  K  in  T.  Ota  et  al.  [167]  model  05  K  and  145 
in  the  first  and  second  models  of  R.  Suwanwarangkul  et  al.  [61  ], 
respectively.  Therefore,  a  strict  temperature  prediction  distri¬ 
bution  within  SOFCs  is  essential  for  anticipating  and  optimizing 
the  cell  performances  more  accurately. 


3.5.3.  Momentum  model 

Momentum  balance  for  the  fuel  in  both  TSOFC  and  PSOFC  yields: 
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Pressure  gradient  in  the  anode  and  cathode  side  is  important  to 
the  electrochemical  activities  of  the  fuel  cell.  It  was  found  that  by 
ignoring  the  anode  pressure  gradient  in  the  model  can  make  and 
error  in  calculating  of  concentration  overpotential  by  about  20% 
at  a  current  density  and  temperature  of  5000 Am-2  and  1073  K, 
respectively  [66]. 


3. 5.3.1.  Different  fuels  fed  SOFCs.  One  of  the  most  brightened  advan¬ 
tages  of  solid  oxide  fuel  cells  over  the  other  types  of  fuel  cells  is  that 
hydrogen  and  carbon  dioxide  are  used  as  fuel  in  the  cell.  This  means 


However;  it  is  not  clear  which  of  these  reactions  predominates 
in  a  SOFC  running  on  H2S. 

Some  researchers  like  Monder  et  al.  [84],  Liu  et  al.  [163],  Slavov 
et  al.  [164]  and  others  in  [84,163,164,180-185]  used  H2S  as  a  fuel 
in  SOFCs. 

Monder  et  al.  [84]  performed  an  isothermal  steady  state  two- 
dimensional  model  of  a  planar  SOFC  fed  by  H2S.  In  this  paper,  the 
performance  of  fuel  cell  fed  by  H2S  is  not  clear  in  comparison  to 
other  type  of  fuels. 

Ahmed  et  al.  [  1 61  ]  were  the  first  group  that  compared  the  effect 
of  different  fuels  on  SOFC  performance.  A  cross-flow  monolithic 
solid  oxide  fuel  cell  honeycomb  (MSOFC)  structure  with  alternat¬ 
ing  layers  of  anode,  electrolyte,  cathode,  and  interconnect  fed  with 
different  fuels,  humidified  hydrogen  and  coal  gas,  and  air  as  oxi¬ 
dants  was  modeled.  Results  of  the  model  with  the  simulated  coal 
gas  as  the  fuel  and  air  as  the  oxidant  showed  that  the  MSOFC  fed 
with  coal  gas  reached  a  same  cell  voltage  as  a  cell  voltage  of  a  hydro¬ 
gen  fed  MSOFC  but  the  coal  gas  MSOFC  was  expected  a  lower  net 
power  and  a  lower  average  current  density. 

Suwanwarangkul  et  al.  [51],  compared  the  performance  of 
TSOFC  using  by  biomass-derived  synthesis  gas  as  a  fuel  with  the 
previous  model  (Suwanwarangkul  et  al.  [61])  for  a  cell  operating 
with  humidified  H2.  The  results  depicted  that  the  cell  performance 
for  operating  with  humidified  H2  is  more  than  that  obtained  with 
the  synthesis  gases  because  higher  average  temperatures  within 
the  cell  were  reached  due  to  the  higher  energy  content  of  humidi¬ 
fied  H2.  On  the  other  hand,  due  to  lower  temperature  distribution 
within  the  cell,  the  cell  average  temperature  fed  with  synthesis  gas 
is  about  100°C  lower  than  that  of  humidified  hydrogen. 

Aloui  et  al.  [69]  presented  an  analytical  modeling  of  polarization 
in  a  solid  oxide  fuel  cell  using  biomass  syngas  product  as  fuel.  The 
performance  of  fuel  cell  fed  by  biomass  syngas  was  compared  with 
the  fuel  cell  fed  by  pure  hydrogen.  It  was  showed  that  a  SOFC  using 
the  direct  oxidation  of  syngas  fuel  (mixture  of  CO  and  H2)  is  more 
powerful  than  that  using  only  H2  fuel. 

Ni  et  al.  [  1 86]  observed  that  the  performance  of  CH4  fed  SOFC  can 
promote  at  high  operating  temperature,  as  all  the  overpotentials 
decrease  with  increasing  temperature.  This  behavior  is  different 
for  H2  fed  SOFC  performance,  because  by  increasing  the  fuel  cell 
temperature  the  concentration  overpotential  also  increases. 

P.  Dokmaingam  et  al.  [77],  performed  a  mathematical  two- 
dimensional  steady  state  model  of  an  indirect  internal  reforming 
tubular  solid  oxide  fuel  cell  fed  by  four  different  primary  fuels,  i.e., 
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methane,  biogas,  methanol  and  ethanol.  The  objective  of  this  paper 
was  to  investigate  the  effect  of  fuel  type  on  the  thermal  coupling 
between  internal  endothermic  reforming  with  exothermic  electro¬ 
chemical  reactions  and  system  performance.  It  was  observed  that 
indirect  internal  reforming  SOFC  fed  by  methanol  has  the  greatest 
performance  in  terms  of  electrical  efficiency  and  temperature  dis¬ 
tribution  along  the  system.  Based  on  this  study,  C02  in  biogas  affects 
the  system  performance  significantly.  The  electrical  efficiency  can 
be  improved  by  removal  of  some  carbon  monoxide  from  biogas; 
however,  a  greater  temperature  gradient  is  expected. 

Brett  et  al.  [187]  were  interested  in  studying  on  methanol  fed 
SOFCs  with  an  intermediate  range  of  temperature  between  500  and 
600  °C.  A  lower  operating  temperature  induces  a  more  rapid  start¬ 
up  and  shut-down,  decreases  the  cell  temperature  fluctuation  and 
then  thermal  stress  in  the  ceramics.  In  addition,  at  a  low  operating 
temperature  the  dimensional  thickness  of  the  electrolyte  can  be 
reduced,  giving  almost  the  same  results  by  promoting  the  ionic  con¬ 
ductivity  of  the  electrolyte  at  lower  operating  temperatures  [159]. 
Moreover,  fueling  an  SOFC  directly  by  methanol  ignores  the  inter¬ 
nal  reforming  process  with  large  thermal  gradients  within  the  fuel 
cell. 

Arpornwichanop  et  al.  [188]  studied  the  performance  of  an 
anode-supported  SOFC  fed  by  ethanol  under  a  direct  internal 
reforming  and  isothermal  condition.  Electrochemical  mode  includ¬ 
ing  all  polarization  losses  was  investigated  in  details.  It  was 
illustrated  that  cell  temperature,  the  ratio  of  steam  to  ethanol  and 
the  inlet  feed  flow  substantially  affect  the  electrical  efficiency  of 
SOFC.  In  this  study,  the  performance  of  fuel  cell  fed  by  ethanol  was 
not  compared  with  other  fuels  fed  solid  oxide  fuel  cell. 


4.  Literature  on  the  SOFC  configuration 

Tubular  and  planer  solid  oxide  fuel  cells  are  discussed  generally 
in  Section  2.  These  two  types  of  SOFCs  are  common  configura¬ 
tions  used  in  the  literatures.  Some  researchers  were  interested  in 
modeling  both  T-SOFC  and  P-SOFC.  Ferguson  et  al.  [50],  studied  a 
three-dimensional  mathematical  model  on  both  tubular  and  plan¬ 
ner  SOFCs  with  considering  a  local  distribution  of  the  electrical 
potential,  temperature  and  concentration  of  the  chemical  species 
inside  the  both  configurations.  By  comparing  both  tubular  and  pla¬ 
nar  solid  oxide  fuel  cells  a  less  ohmic  loss  and  higher  efficiency  for 
the  former  geometry  was  observed. 

Stiller  et  al.  [28]  like  Ferguson’s  group  [50]  were  interested  in 
both  planar  and  tubular  solid  oxide  fuel  cells  modeling.  However; 
the  study  was  focused  more  on  effects  of  parameters  on  both  fuel 
cell  configurations  with  a  2D  steady-state  planar  and  tubular  SOFCs 
fed  by  methane.  The  model  was  involved  of  ohmic  resistance,  con¬ 
vective,  conductive  and  radiation  heat  transfers,  chemical  reactions 
and  mass  conservations.  The  required  boundary  conditions  and 
input  parameters  of  the  models  were  discussed.  The  objective  of 
the  study  was  to  observe  the  impact  of  inlet  air  flow  temperature, 
pressure  ratio,  air  flow  rate  and  anode  gas  recycling  on  both  con¬ 
figurations.  It  was  illustrated  that  the  cycle  balance  is  different  in 
both  systems,  as  the  tubular  fuel  cell  required  a  lower  air  inlet 
temperature.  Both  fuel  cell  systems  achieved  above  65%  electric 
efficiency. 

A  complete  comparison  between  tubular  and  planar  solid  oxide 
fuel  cell  configurations  was  presented  by  Zhang  et  al.  [156].  A 
three-dimensional  steady-state  numerical  model  for  both  config¬ 
urations  was  described.  This  paper  conducted  detailed  analyses  of 
the  electric  polarizations  for  both  configurations  and  presented  the 
characteristics  of  the  potential  distribution  in  the  electric  com¬ 
ponents  for  the  planar  co-flow  SOFC.  The  results  were  in  a  good 
agreement  with  the  experimental  data.  Based  on  the  results,  for 
planar  SOFC,  polarizations  results  in  the  electrolyte  ohmic  loss  and 


the  activation  polarization  but  for  the  tubular  SOFC,  the  cathode 
ohmic  loss  has  the  greatest  contribution  to  the  sum  of  polariza¬ 
tions  while  the  activation  polarization  is  the  next  highest  loss.  It 
was  also  shown  that  the  ohmic  heat  predominately  depends  on  the 
current  density  in  the  planar  SOFC  while  in  tubular  configuration 
the  ohmic  heat  is  produced  from  the  cathode.  The  relative  ohmic 
heat  contribution  to  the  total  heat  is  higher  for  the  tubular  than 
for  the  planar  configuration,  and  for  the  planar  SOFC,  this  is  more 
pronounced  in  the  counter-flow  design  than  in  the  co-flow  design. 
Before  Zhang  et  al.  group  [156]  the  effect  of  different  flow  manifold¬ 
ing,  i.e.,  cross-co,  or  counter-flow  were  considered  by  Achenbach 
[  1 1 6].  It  was  showed  that  the  highest  fuel  cell  efficiency  is  obtained 
for  counter-flow,  the  most  uniform  current  density  distribution  is 
reached  for  co-flow  and  the  largest  temperature  gradients  in  the 
solid  structure  occur  for  cross-flow.  Ferguson  et  al.  [50]  compared 
the  counter  flow  and  co-flow  for  planar  geometry  in  their  model 
and  revealed  that  the  counter  flow  for  hydrogen  as  a  fuel  is  the  most 
suitable  flow  configuration.  A  comparison  between  counter-flow 
and  co-flow  was  performed  by  Aguiar  et  al.  [  1 08  ].  It  was  presented 
that  for  the  counter-flow  configuration  the  temperature  increases 
along  the  air  flow  direction  with  the  maximum  at  the  outlet  and 
for  the  co-flow  configuration  the  temperature  profile  reaches  to 
its  maximum  value  near  the  fuel  flow  outlet.  In  this  system  co¬ 
flow  configuration  is  the  more  efficient  one  due  to  the  high  air 
temperature. 


5.  Novel  SOFC  models 

As  a  kind  of  nonlinear  multi-input-multi-output  (MIMO)  sys¬ 
tem,  SOFC  is  hard  to  model  by  the  traditional  methodologies. 
Therefore;  some  researchers  have  attempted  to  establish  novel 
fuel  cell  models  by  statistical  data-driven  approach.  Data-driven 
approach  is  based  on  analyzing  the  data  about  a  system,  in  partic¬ 
ular  finding  relationships  between  the  input  and  output  variables 
of  the  system  without  explicit  knowledge  of  its  physical  behavior. 
The  artificial  neural  network  model  is  a  novel  model  which  is  very 
useful  for  rapid  calculation  of  SOFC  e.g.  dynamic,  static  and  control 
simulations.  The  Hammerstein  models  is  also  a  novel  method  for 
modeling  SOFCs  which  include  a  nonlinear  block  followed  by  the 
linear  system  and  are  useful  to  model  static  and  dynamic  systems 
of  the  solid  oxide  fuel  cell. 

Some  researchers  have  attempted  to  establish  novel  SOFC  mod¬ 
els  [189,190,194-196]. 

Jurado  [149]  developed  a  multivariable  Hammerstein  model  of 
a  SOFC,  which  is  applicable  for  small  signal  and  dynamic  stabil¬ 
ity  studies.  The  Hammerstein  models  include  a  nonlinear  block 
followed  by  the  linear  system  and  are  useful  to  model  static  and 
dynamic  systems  of  the  solid  oxide  fuel  cell.  Dynamic  behavior 
of  temperature  and  output  voltage  losses  was  explained  by  con¬ 
sidering  the  electrochemical  and  thermal  equations  in  this  paper. 
Non-iterative  algorithms  for  the  identification  of  multivariable 
Hammerstein  system  were  studied.  The  main  debate  in  the  deriva¬ 
tion  of  the  results  was  the  representation  of  the  system  using  base 
functions. 

Bo  Huo  et  al.  [191  ]  reported  a  dynamic  modeling  of  SOFC  using 
Hammerstein  Model  for  dynamic  simulation  and  control.  The  static 
nonlinear  of  the  Hammerstein  model  was  modeled  by  a  radial  basis 
function  neural  network  (RBFNN),  and  the  linear  part  was  mod¬ 
eled  by  an  autoregressive  with  exogenous  input  (ARX)  model.  In 
this  study,  a  new  gradient  descent  algorithm  was  used  to  antici¬ 
pate  the  hidden  centers,  the  radial  basis  function  widths  and  the 
connection  weights  of  the  RBFNN.  For  anticipating  the  parameters 
and  the  orders  of  the  ARX  model,  the  least  squares  (LS)  algorithm 
and  Akaike  Information  Criteria  (AIC)  were  derived,  respectively.  By 
comparing  the  Hammerstein  model  and  the  RBFNN  model  exper- 
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Model  including  conduction/convection/radiation/electrochemical  heating/ 
diffusional  heating  (T.  Ota  et  al.  [54]) 

Model  with  radiation  (R.  Suwanwarangkul.  et  al.  [42]) 

Model  without  radiation  (R.  Suwanwarangkul.  et  [42]) 


Fig.  6.  Comparison  of  cell  temperature  distribution  between  three  different  models. 

imentally,  a  better  performance  for  the  Hammerstein  model  was 
reported  in  this  study. 

Wu  et  al.  [192]  presented  a  nonlinear  offline  model  of  SOFC 
used  a  radial  basis  function  neural  network  based  on  genetic  algo¬ 
rithm.  The  objective  of  the  work  was  to  optimize  the  parameters 
of  radial  basis  function  neural  networks.  The  model  was  validated 
by  the  simulations.  The  simulation  results  illustrated  that  radial 
basis  function  neural  network  based  on  genetic  algorithm  method 
is  superior  to  the  conventional  back  progration  neural  network 
in  predicting  the  stack  voltage  with  different  temperature  basis 
parameter  neural  network  in  anticipating  the  fuel  cell  voltage 
in  case  of  different  temperature.  Although  many  other  operating 
parameters  are  also  effective  on  fuel  cell  performance,  only  current 
density  and  temperature  was  considered  in  this  model.  Wu  group 
[193]  also  studied  the  different  flows  effect  on  the  performance 
of  SOFC  by  developing  an  adaptive  neural-fuzzy  inference  system 
(ANFIS).  In  order  to  identify  linear  and  nonlinear  parameters  in  the 
ANFIS,  a  hybrid  learning  algorithm  combining  back  propagation 
and  least  squares  estimate  was  adopted.  The  results  were  com¬ 
pared  by  simulations  and  the  simulation  results  showed  that  the 
ANFIS  model  can  efficiency  predict  the  dynamic  behavior  of  SOFC. 
However,  the  model  was  not  included  the  temperature  dynamic 
response  and  the  SOFC  fuel  processor  was  also  not  considered. 

Milewski  et  al.  [190]  found  that  the  Artificial  Neural  Network  can 
be  used  to  simulate  the  fuel  cell  behavior  only  by  utilizing  available 
experimental  data.  The  cell  voltage  was  predicted  by  the  inputs  of 
the  model  which  were  current  density,  temperature,  fuel  volume 
flow  density,  and  oxidant  volume  flow  density.  The  results  reveal 
that  the  ANN  can  be  efficiently  used  for  the  single  SOFC  with  rela¬ 
tively  high  accuracy.  Table  7  shows  some  researchers  that  focused 
on  novel  models. 

6.  Conclusion 

Mathematical  modeling  is  an  essential  tool  in  design  of  fuel  cell 
systems,  as  it  is  important  to  understand  the  response  of  a  cell  stack 
under  normal  or  transient  conditions.  Such  models  can  provide  a 
picture  of  cell  stress,  potential,  current  density,  and  temperature  as 
functions  of  position  and  time  for  various  cell  configurations  and 
operating  conditions,  and  be  used  to  examine  the  effects  of  changes 
on  one  or  more  variables  and  the  relative  system  sensitivity  to  rel¬ 
evant  design  parameters.  A  model  is  also  useful  in  predicting  the 
effects  of  altering  process  variables  and  in  using  that  information  to 


optimize  cell  performance.  Many  research  studies  have  been  done 
in  this  regard.  However,  still  many  aspects  of  fuel  cell  modeling  are 
remained  for  studying  (Fig.  6). 

In  this  study,  mathematical  modeling  of  SOFCs  for  both  tubu¬ 
lar  and  planar  configurations  was  reviewed  into  five  subsystems 
considering  the  factors  like  polarization  losses,  mass/energy/ 
momentum  conservations,  diffusion  through  porous  media,  and 
electrochemical  phenomena  in  the  PEN  region  and  shift/reforming 
reactions  inside  them  which  is  novel  in  this  regard.  Also,  using  vari¬ 
ety  of  fuels  fed  to  the  SOFCs  was  discussed  and  their  effect  on  the 
system  was  compared  briefly.  A  short  review  of  solid  oxide  fuel  cell 
configurations  and  different  flow  manifolding  were  also  presented 
in  this  review  study.  For  the  first  time  in  this  area,  novel  models 
based  on  statistical  data-driven  approach  existing  in  the  literatures 
were  studied. 

In  this  section  we  point  to  the  subjects  that  need  further 
research. 

Fuel  cell  temperature  and  voltage  strikingly  depend  on  electro¬ 
chemical  reaction.  The  electrochemical  oxidation  of  both  hydrogen 
and  carbon  monoxide  occurs  on  the  triple  phase  boundary,  whereas 
the  proportion  of  these  two  reactions  is  still  unknown.  CO  oxidation 
influences  the  overall  performance  of  the  cell  and  on  the  other  hand, 
polarization  resistances  impress  voltage  directly.  In  many  studies 
CO  oxidation  at  electrochemical  reaction  is  neglected  or  polariza¬ 
tion  losses  are  not  completely  accounted.  However,  it  is  not  still 
quite  clear  that  by  making  some  simplified  assumptions  for  the 
electrochemical  modeling  how  they  can  influence  the  accuracy  of 
the  model.  What  the  difference  between  observed  experimental 
values  and  the  results  of  the  mathematical  modeling  is.  Calculating 
the  errors  even  for  small  changes  means  that  there  is  something 
we  do  not  understand  and  need  to  account  for  it.  Therefore,  such 
phenomena  should  be  deeply  studied  in  both  mathematical  mod¬ 
eling  and  experimentation  fields  to  perform  a  more  accurate  and 
validated  model. 

Very  limited  studies  focused  on  how  the  SOFC  microstructure 
influence  the  chemical  reactions  and  mass  transfer  in  the  porous 
SOFC  electrode.  Generally,  length  of  TPB,  porosity,  particle  size  and 
thicknesses  of  electrodes  and  electrolytes  has  conflicting  effect  on 
the  PEN  performance.  Thus,  optimal  design  for  PEN  in  order  to 
improve  the  PEN  efficiency  is  significant. 

Diffusion  through  porous  media  significantly  affects  the  perfor¬ 
mance  of  solid  oxide  fuel  cells.  Fuel  cell  voltage  and  polarization 
losses  depend  on  the  concentration  of  components  inside  the  triple 
phase  boundary.  In  addition,  diffusional  heat  transfer  influences 
the  fuel  cell  temperature.  Therefore,  it  is  necessary  to  consider 
diffusion  in  models  in  order  to  have  an  accurate  prediction  of 
fuel  cell  behaviors.  Although  several  research  has  been  done  for 
presenting  a  beneficial  diffusional  model  (FM,  SMM,  DGM)  at  differ¬ 
ent  conditions,  still  more  investigation  requires  to  compare  these 
three  models  together  to  minimize  the  error  between  mathemati¬ 
cal  modeling  and  experimental  results. 

Among  the  papers  which  were  studied  in  this  work,  was  shown 
that  some  of  the  models  ignored  radiative  heat  transfer  between 
free  walls  and  especial  between  gas  components.  Many  works 
accounted  radiation  heat  transfer  in  their  model  but  the  knowledge 
of  thermal  radiation  property  for  both  bulk  and  surface  proper¬ 
ties  at  high  temperatures  relevant  to  SOFC  is  still  incomplete  and 
requires  more  study  in  this  regard. 

Electrochemical  heating  is  also  still  not  fully  available  in  open 
literature  and  is  imperative  to  develop  in  the  models.  Ohmic  heat¬ 
ing  must  be  tacked  into  account  in  SOFC  models  due  to  the  large 
ionic  resistance  of  YSZ  compared  to  heat  of  electrolyte  materials 
used  in  other  types  of  fuel  cells.  Moreover,  conduction  heat  transfer 
in  many  papers  was  neglected  while  accurate  determination  of  the 
temperature  distribution  including  the  radiation,  conduction,  diffu¬ 
sional  heat  transfer  and  electrochemical  heating  is  required  since 
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the  material  properties,  chemical  kinetics,  and  transport  proper¬ 
ties  of  materials  used  in  single  cell  and  stack  level  depends  on  the 
temperature. 

Neglecting  the  momentum  balance  by  several  researchers  was 
not  fully  justified.  More  research  studies  require  showing  the  effect 
of  momentum  balance  especially  inside  the  porous  electrodes  on 
solid  oxide  fuel  cell  performance. 

The  pressure  gradient  in  the  electrodes  of  SOFCs  is  simply 
ignored  without  any  justification  in  modeling  studies.  Therefore 
the  effect  of  pressure  distribution  on  mass  transfer  and  thus  on 
overpotential  should  be  more  studied. 

A  major  advantage  of  solid  oxide  fuel  cells  over  other  types 
of  fuel  cells  is  that  a  variety  of  hydrocarbon-based  gases  or  their 
synthesis  derivatives,  such  as  natural  gas,  biomass  and  coal  can 
potentially  be  used  as  fuel  sources.  It  was  shown  that  a  SOFC  using 
the  direct  oxidation  of  syngas  fuel  (mixture  of  CO  and  H2)  is  more 
powerful  than  using  pure  H2  fuel,  and  also  methanol  in  compar¬ 
ison  to  other  fuels  like  methane,  biogas,  and  ethanol  as  fuels,  has 
the  greatest  performance  in  terms  of  electrical  efficiency.  Flowever, 
still  a  complete  research  study  needs  to  present  a  suitable  fuel  for 
promoting  the  SOFCs  performance  and  their  power  efficiency.  In 
addition,  the  effect  of  each  fuel  on  the  important  factors  of  SOFCs 
such  as  cell  temperature,  voltage,  reactions,  and  cell  performance 
is  not  completely  studied. 

In  many  papers  the  effect  of  reformers  and  different  geome¬ 
tries  and  configurations  is  not  clear.  It  is  imperative  to  show  the 
importance  of  the  mentioned  factors  in  the  models. 

The  intelligent  novel  models  avoid  the  complexity  of  the  math¬ 
ematical  model  and  provide  the  input-output  relationship  much 
faster  and  easier.  Based  on  these  models,  valid  control  strategy 
studies  like  predictive  control  and  robust  control  can  be  developed. 
Therefore,  more  research  studies  requires  for  developing  the  novel 
models. 
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